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Abstract 

An effective theory based on ensembles of either regular gauge instantons or 
merons is shown to produce confinement in SU(2) Yang-Mills theory. When the 
scale is set by the string tension, the action density, topological susceptibility and 
low-lying glueball spectrum are similar to those arising in lattice QCD. The phys- 
ical mechanism producing confinement is explained, and a number of analytical 
insights into the effective theory are presented. 
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1 Introduction 



Understanding the mechanism of confinement in quantum chromodynamics (QCD) is a funda- 
mental but elusive challenge. Lattice calculations convincingly demonstrate that QCD yields 
confinement, but do not yet elucidate its mechanism. Although the strong coupling expansion 
produces an area law by disorder already in lowest order, this disorder alone is not the whole 
story as evidenced by the fact that the strong coupling expansion also erroneously yields con- 
finement in U(l) gauge theory. Hence, in this work, we develop a confining effective theory 
based on degrees of freedom arising naturally from QCD. We explore in detail how confine- 
ment arises, and the degree to which this effective theory approximates the results of SU(2) 
lattice gauge theory. 

There has been a long and fruitful history of seeking to understand nonperturbative QCD an- 
alytically by expanding the path integral for the partition function around stationary classical 
solutions and evaluating the fluctuations around these solutions. Central to this approach was 
the discovery of instantons HI that satisfy the classical Euclidean Yang-Mills equations and 
implement the tunneling required between different winding number sectors of the vacuum. 
Since instantons can be transformed to the singular gauge in which the gauge field falls off like 
~ 1/x 3 , a dilute gas of these instantons is an approximate solution to the field equations. The 
contribution of this dilute instanton gas to the potential between two static color charges can 
be calculated analytically, and was shown not to produce confinement [2J. Whereas an instan- 
ton that contributes to two opposite Wilson lines in a Wilson loop does contribute to a linear 
term to the potential, when the integral over color orientations of independent instantons that 
affect only one or the other Wilson line is performed, these uncorrelated contributions do not 
produce a linear term. Hence at separation larger than the instanton size, the potential becomes 
a constant and fails to confine. 

When instantons failed to produce confinement, there was then hope that merons [3], solu- 
tions to the Yang-Mills equations with topological charge 1/2, would do so. Although a free 
meron has an action that diverges logarithmically with the volume, a pair of merons has a fi- 
nite action that depends logarithmically on the separation. If the coupling is strong enough that 
the logarithmic potential is weaker than the entropy associated with the Euclidean separation 
(~ r 3 = e 31nr ), meron pairs could dissociate and form a meron gas. Since meron fields only 
fall off as ~ 1/x, they are expected to be more effective in disordering than singular gauge 
instantons and thus potentially be a more effective mechanism for confinement flU[51[6]|. Unfor- 
tunately the same ~ 1/x behavior that makes them a candidate for confinement renders them 
intractable analytically, and as a result they have never been shown by analytic arguments to be 
responsible for confinement. 

Motivated by the quark zero modes associated with instantons and the 't Hooft interaction gov- 
erning light quarks, highly successful instanton liquid models, which had to be implemented 
numerically, were developed that provide clear insight into the role instantons play in chi- 
ral symmetry breaking and reproduce many of the salient features of QCD with light quarks 
B71 [H |H QU. For our present purposes, it is best to think of an instanton liquid as an effec- 
tive theory in which one parameterizes a set of gauge configurations in terms of the positions 
and color orientations of an ensemble of instantons to represent the salient degrees of freedom 
for gluons. Because the field equations are nonlinear, the superposition of a set of instanton 
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solutions is not a classical solution to the field equations. Hence, rather than summing over a 
set of classical solutions, one is approximating the original functional integral for the partition 
function over gauge fields by performing the integral over these effective degrees of freedom, 
for example by using the Metropolis algorithm to sample configurations. Physically, this sam- 
pling will favor the appropriate configurations that are energetically and entropically optimal. 
Because the gauge fields of singular gauge instantons fall off as ~ 1/a; 3 , a random instanton 
liquid in which the positions and color orientations of instantons are chosen randomly does not 
differ qualitatively from one generated by more sophisticated Metropolis sampling, so most of 
the qualitative phenomena occur in a random instanton liquid. As in the case of the dilute in- 
stanton gas, Wilson loops in a random instanton liquid also show no confinement at distances 
larger than the instanton size. 

In the present work, we consider a more general class of effective theories based on ensembles 
of regular gauge instantons or merons, whose fields fall off like ~ l/x, so that the resulting 
superposition of gauge fields is not close to a classical solution. Throughout this work, we will 
refer to regular gauge instantons or merons as pseudoparticles. The term instanton will always 
mean a regular gauge instanton, and if we wish to refer to a singular gauge instanton, we will 
always do so explicitly. Although it is analytically intractable to perform the resulting integral of 
the partition function over the pseudoparticle collective coordinates, it is straightforward to do 
so numerically using the Metropolis algorithm. Whereas most random choices of pseudoparticle 
fields create a large background field yielding an unphysically high action density, Metropolis 
sampling selects the particular configurations in which the collective coordinates are correlated 
in such a way as to produce very low action density. In this way, one can achieve the balance 
of energy and entropy in ensembles based on regular gauge pseudoparticles that is conceptually 
equivalent to what has been done in the past with singular gauge instantons. 

This more general class of theories based on pseudoparticles combines both confinement and 
the successes of previous singular gauge instanton liquid models. As will be shown below, color 
correlations are required to superpose the long range pseudoparticle gauge fields to create the 
small background fields whose low action allows them to dominate the path integral. These also 
produce the correlations necessary to yield an area law for large Wilson loops and thus produce 
confinement. The short distance "spikes" in the gluon field arising from the short distance be- 
havior of small merons or instantons behave just like the spikes in conventional singular gauge 
instanton liquid models, and thus produce the physics associated with the 't Hooft interaction 
and chiral symmetry breaking. Pseudoparticle ensembles were first used by the present authors 
to study meron ensembles IfTTTl . where there is no singular gauge alternative. However, in the 
context we have just described, it is clear that one should also consider regular gauge instantons. 
Indeed, since the gauge for each instanton is defined relative to its position, the fields obtained 
by superposing N singular gauge instantons is drastically different from the one obtained by 
transforming each of these instantons to the regular gauge and superposing the resulting fields. 
Neither is a classical solution of the Yang-Mills equations, and neither is preferred on physical 
grounds. Hence, in exploring the full physics based on classical solutions, it is essential to study 
both alternatives. Although we will not develop it further in this work, Ref. IfTTTl discusses how 
one may expect center symmetry to be realized in pseudoparticle effective theories. A closely 
related effective theory lfT2l superposes pseudoparticle fields with dynamically determined am- 
plitudes and generalizes to finite temperature, where a clear crossover in the Polyakov loop is 
observed from a nearly vanishing value at zero temperature to a magnitude close to unity at high 
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temperature. 



The primary goals of this work are to understand the mechanism by which this effective the- 
ory produces confinement and other observed properties of SU(2) gauge theory, and explore its 
success and limitations in quantitatively approximating the results of lattice QCD calculations 
of the action density, topological susceptibility, and glueball masses. Although we will discuss 
both instanton and meron pseudoparticles in detail and give extensive examples for both, the 
primary focus in comparing with the lattice QCD spectrum is on the confining instanton theory. 
We will show that the confining instanton theory has essentially all the advantages of previous 
random and interacting singular gauge instanton models, with the added feature of incorporat- 
ing confinement. To obtain insight into the effective pseudoparticle theory, we will strongly 
emphasize analytic arguments wherever possible, as indicated in the outline below. 

The outline of this paper is as follows. Section 2 describes the details of the pseudoparticle 
ensembles and some relevant analytical properties. In particular, since one might be concerned 
that the long range fields of pseudoparticles are necessarily large, we present a simple example 
of how sets of merons can be chosen such that the leading order 1/x 4 term in the action at 
large distances vanishes identically. Section 3 describes random and correlated pseudoparticle 
ensembles. For pedagogical purposes, we begin with an ensemble in which the pseudoparti- 
cle positions and color orientations are selected randomly without correlations. Although this 
ensemble has unphysical features, it enables us to make instructive analytical arguments. We 
then allow Metropolis dynamics to produce the appropriate color and spatial correlations. The 
resulting gauge fields are shown to have a proper thermodynamic limit and to be composed of 
small background fields plus "spikes" at the locations of pseudoparticles. By separately apply- 
ing Metropolis dynamics to only color parameters and then both color and spatial parameters, 
we show that the dominant correlations are in the color degrees of freedom. Section 4 addresses 
Wilson loops and confinement. When scaled appropriately, Wilson loops are shown to lie on a 
universal scaling curve, which is nearly independent of the pseudoparticle size and the coupling 
constant and which shows a clear area law indicating confinement. The behavior of Wilson 
loops in both the small size and large size limits is explained analytically. We also show that the 
distribution of Wilson loops follows the diffusion behavior of Ref. [fT3l and that Wilson loops 
in higher representations exhibit Casimir scaling, both of which are in agreement with lattice 
QCD. Section 5 presents calculations of physical observables. Since the theory is confining, 
we can set the scale by the measured string tension and compare unambiguously with lattice 
calculations in which observables are also specified in terms of the measured string tension. 
The two most robust observables are the gluon condensate and the topological susceptibility, 
and both instanton and meron ensembles are shown to yield values in units of the string ten- 
sion in qualitative agreement with lattice QCD. We show that for many operators of interest, 
the effective theory produces Euclidean correlation functions with the proper physical behavior. 
The correlation function of the topological charge density displays the proper negativity behav- 
ior, and operators constructed from field-strengths and Wilson loops enable us to measure the 
masses of + , l - , 1 + , 2 + , and 2~ glueball states. Aside from an overall scale factor which we 
discuss, this glueball spectrum is in qualitative agreement with lattice QCD. The conclusions 
are summarized in the final Section 6. 
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2 Ensembles of Pseudoparticles 



2.1 Pseudoparticle Properties and Definition of Pseudoparticle Ensem- 
bles 



In this work, we explore the idea that pseudoparticles, i.e., merons or regular gauge instantons, 
are the essential degrees of freedom in SU(2) Yang-Mills gauge theory by writing the partition 
function as a path integral of an effective action depending on the positions and color orienta- 
tions of an ensemble of pseudoparticles, 



Z 



^S[A(zi,hi)] 



We identify the effective action S with the standard (Euclidean) action 
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associated with the fields in the space-time region V. Dual field-strength and topological charge 
density are defined by 
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The gauge field for a pseudoparticle in Lorentz gauge with its center at the origin, after appro- 
priate choice of the coordinate system in color space and after regularization of the singularity, 
is given by 

Va^tu-^u & 



a^x) 
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with £ = 1 denoting merons and £ = 2 instantons. Color and space-time dependence are 
correlated via the 't Hooft tensor 
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For instantons and for merons of vanishing size p, a^(x) is a solution of the Euclidean classical 
field equations HI El. The gauge fields of anti-merons or anti-instantons differ in sign if one of 
the space-time indices is 0. The field-strength of the pseudoparticles, its dual, the action density, 
the topological charge density and the topological charge are given by 
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U= ^J d 4 ^ (x) = ±i. (9) 

Pseudoparticles and anti-pseudoparticles differ in the sign of the topological charge and its den- 
sity. We note the difference in the asymptotic behavior of the field- strength of instantons and 
merons. The meron field- strength decays asymptotically as l/x 2 , giving rise to an infrared 
logarithmic singularity in the action. For vanishing size, the action of a meron is also logarith- 
mically divergent in the ultraviolet, so for merons, the size parameter p acts as an ultraviolet 
regulator. 

An important quantity in our study of confinement will be the Wilson loop 

W = ^trl^Pexpig (j) A^(x)dx^\ . (10) 



The integral in (flOl) is ordered along the closed path C. Here we compile the relevant result 
for Wilson loops in the field of single pseudoparticles. Pseudoparticle fields [Eq. ©] in a 
"spacelike" plane with x 4 = are of the form 

£ x x cr 
a — - 



2 x 2 + p 2 



Analytical expressions can be derived if the center of the pseudoparticle fields is in the plane of 
the Wilson loop. In this case, for a loop in the (x±, x 2 ) plane, 

_ £ x 2 £ x x 

2 x{ + x\ + p 2 2 x\ + x\ + p 2 

and the field is abelian pointing in the 3 direction in color space. For calculation of the Wilson 
loop we can use Stokes' theorem with 



(rota 3 ) 3 = ~ 



i p 2 



2 (x 2 + x 2 + p 2 ) 2 

and obtain for a circular Wilson loop (radius r) with the center of the circle and of the pseu- 
doparticle coinciding 

f r 2 
a M dx M = / dcrrota = £7r — -a 3 . 

J-e r + P 

In the limit of vanishing size, 

rota 3 = £7r5^(x)<X3, p^O 

and therefore the Wilson loop in a plane containing the singularity of the pseudoparticle field 
(position x) becomes 

W(x) = (-l) 59(x) , 0(x) = 1 (0) if x is inside (outside) the loop. (11) 

Instantons generate a trivial Wilson loop. This is true everywhere, not only if the center is 
located on the plane of the Wilson loop. For p = the instanton field is a pure gauge, 

= ifij^nj (12) 
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with 



x 4 + icrx 



fi/ = - 2 , (13) 

and therefore loops not passing through the singularity generate 

W lnst EE 1. 

The Wilson loop of merons cannot be evaluated in closed form if the plane of the loop does not 
contain the singularity. 

The ensembles to be considered in this study contain field configurations obtained by superpo- 
sition of pseudoparticles and anti-pseudoparticles of fixed and equal number N P /2, 

N P 

A^x) = J2 M*K(z - z{i))hT\i). (14) 
i=i 

Such a configuration is specified by the position of the centers z(i) and their color orientations 

h(i) = ho(i) + ih(i) • a , h 2 Q (i) + h 2 (i) = 1. (15) 

In the ensembles to be discussed, the location of the merons is restricted to a hypercube 

-1 V = 16. 

Our standard choice for the meron size and coupling constant is 

p = 0.16, g 2 = S2. (16) 
We will see later that after rescaling, our results have very little dependence on these values. 

2.2 Finite Action Density of Pseudoparticles 

The infrared divergence in the action of a meron may appear as an obstacle in their use as 
building blocks of field configurations with (infrared) finite action. A similar concern arises in 
the superposition of (regular gauge) instantons. However, here we show that with an appropriate 
choice of the color orientation, the action of a cluster of merons of finite size can be made finite. 
Similar arguments apply for clusters of regular gauge instantons. 

By dimensional arguments, the leading term of the action density of a collection of merons 
behaves asymptotically as 

and therefore leads to a total action that is logarithmically divergent in the infrared. To calculate 
the coefficient S^y we proceed by averaging over the orientation. Useful identities are: 



1 / .... 1 
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We introduce the notation 



Q 1 (]ni -^ 4 



Vatiu if 2-th pseudoparticle is a meron , 
fja/xv if i-th pseudoparticle is an anti-meron . 



and find 

Veflp 



Vepp 

&ik&jl&ee'&ff' + 8il8jk8ef'8fe' + (1 — ^ijW&ij&kl&effie 1 f 



After a tedious calculation, we arrive at the following expression for the coefficient S^: 
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(18) 



The color orientation of 4 merons can be chosen such as to make the leading asymptotic contri- 
bution to the action vanish. This corresponds to a neutral object. From the trivial identity 
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and by comparison with the representation of the color matrices u 

-6 bc (h 2 -h 2 )+ 2h h a e abc + 2h b h c 



u 



(19) 



we read off that 4 merons with the color orientation of the i-th meron given by hi = 1 yield a 
vanishing l/x 4 term in the asymptotic expansion of s. Such configurations exhibit confinement 
of their building blocks, analogous to the role of the gluons, since a complete dissociation of 
such a pseudoparticle cluster is accompanied by an infinite increase in the action. 

The pseudoparticle ensembles have been generated by Monte Carlo sampling of the action in 
the path integral, Eq. (QQ). In each step of a Metropolis update, the position and color orientation 
of a given meron are tentatively changed, the induced changes in the action density are evaluated 
at a set of mesh points distributed over the whole volume, and the configuration is accepted or 
rejected based on the global change in action. The long range nature of the meron fields makes 
the changes extend throughout the whole system. 
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3 Stochastic and Correlated Pseudoparticle Ensembles 



Our goal is to calculate vacuum expectation values of observables via the path integral over the 
collective coordinates of our effective theory with pseudoparticle degrees of freedom 

(0) = ^J dzd^e-^^OlAiz,, hi)], (20) 

where Z is the partition function defined in Eq. (OQ). As discussed in the introduction, it is 
useful to proceed in steps. We begin with an ensemble in which the pseudoparticle positions 
and color orientations are selected according to a uniform distribution. This corresponds to 
the strong coupling limit, in which g 2 is infinite so that the action does not affect the weight 
of each configuration. It will be fruitful to think of the present effective theory in the context 
of corresponding lattice calculations, and this stochastic ensemble will be comparable to the 
strong coupling limit of lattice QCD. As in lattice QCD, various quantities can be computed 
analytically in the strong coupling limit. The comparison with the stochastic ensembles will 
be a cornerstone for understanding the significantly more complex structure of the dynamically 
correlated ensembles. 

Dynamical correlations are included using the Metropolis algorithm to sample variables with a 
probability distribution given by the action. Thus, in the usual way, to calculate the path integral 
J dxe~ s ^O(x), we change x with a microreversible change Ax, calculate the corresponding 
change in action, AS, and accept the new value x + Ax with probability P = min (1, e~ A5 ). 

The crucial step beyond the stochastic ensemble is to include dynamical correlations in the 
pseudoparticle color orientations. Update of pseudoparticle positions is much less important, 
but will always be included unless we want to highlight the difference. Thus both the positions 
z(i) and the parameters h(i) specifying the color orientations of all the pseudoparticles are equi- 
librated using Metropolis updates. We will show that these dynamical color correlations make a 
qualitative change in the physical observables. In order to allow for a variable ultraviolet scale, 
we will also consider ensembles of pseudoparticles with a uniformly distributed pseudoparticle 
size p. 



3.1 Action Density 

In the stochastic ensemble, the action density can be evaluated analytically by assuming that the 
average over the color orientations satisfies 

u ac {i)u bc ( 3 )^5, l3 5 ah . 

This leads to the following expression: 

N P 



As re 



W) W) + E [(4(0 4(0) K(j) Ku)) - (Al{i) A b v {ij) (Ai{j) Ai{j))] 



Using standard identities for the 't Hooft symbol, we obtain, cf. Eq. © 

g 2 Aft) A%i) = - ? - p2]2 (<W (x - z(z)) 2 - (x - (x - z u (i)) ) 



which in the center of the cube (x = 0) yields 
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(21) 



Furthermore, with the action density of a single meron or instanton, cf. Eq. ©, 

6£ 2 
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we derive for the action density of the gauge field generated by N P pseudoparticles 
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(23) 



Performing the ensemble average we obtain for a system of N M merons 

g 2 (s(x = 0)> = 4.5 N M + 1.03 N M (N M - 1). 



(24) 



Fig.[T]shows the action density for the stochastic and for dynamically correlated ensembles. We 
observe a qualitatively different behavior of the action density as function of the position in the 
cube for stochastic and dynamically correlated ensembles. In agreement with the expression 
(123b . the action density of the stochastic ensemble peaks in the center of the cube and decreases 
by a factor of 2 when approaching the surface of the cube. The dynamically correlated en- 
sembles, on the other hand, are approximately constant for \x^\ < 0.8 and thus seem to respect 
translational symmetry. This central region is used for the evaluation of observables. The values 
of the action density of stochastic and dynamically correlated ensembles differ by two orders of 
magnitude. The value of the action density of the stochastic ensemble in the center of the cube, 

g 2 (s(x = 0)) = 2.6 • 10 5 , 



is, within 1%, correctly predicted by Eq. (|24]>. 

The one-pseudoparticle contribution to the action density in (124)) arising from the abelian part 
of the field-strength yields an intensive contribution to the action density in the thermodynamic 
limit, up to a logarithmic correction. In contrast, in the thermodynamic limit, assuming a fi- 
nite density of pseudoparticles, the non-abelian contributions to the field-strength are dominant 
and give rise to an action density increasing linearly with the volume for sufficiently small 
p. Stochastic cancellations are not sufficient to achieve volume independence. For a cube of 
side-length 2L with L — » oo, the action density of N M merons at the center of the cube is 



(s(x = 0)) w 1.29 



N 



M 



g 2 L*- 



(25) 



In our calculation, all interference terms have disappeared and the same result applies for an 
arbitrary mixture of a total of N M merons and anti-merons; furthermore, the same number of 
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Figure 1: Action density g 2 (s(x)) along the x-axis for the stochastic (rescaled) (green x) and 
dynamically correlated ensembles with stochastic (blue *) and with dynamical positions (red 
+). Dynamical correlations decrease the value at the center of the cube from 261 000 to 1700. 

regular instantons yields a 16 times larger action density. Clearly, the expression (1251) is not 
compatible with a proper thermodynamic limit. 

For comparison, we refer to the action density in the standard dilute instanton liquid. In this 
model the building blocks are instantons in singular gauge with the gauge field decreasing 
asymptotically like 

A a ~ — . 

Thus for sufficiently low densities of instantons, nj <C p -4 , the action density is given by 

g 2 (s) = 8n 2 ni (26) 
and yields for Nj instantons in the volume 2 4 

g 2 {s) «4.94JV>. (27) 

The action associated with a field generated by stochastic superposition of 500 merons (regular 
instantons) is about 100 (16 000) times larger than the action of 500 singular instantons. As 
the expression (l24l) shows, the large value of the action density is not generated by pseudoparti- 
cles in the neighborhood, rather it results from the superposition of gauge fields with their l/x 
asymptotic behavior generated by distant pseudoparticles. The contour plot of the action den- 
sity of a single configuration of the stochastic meron ensemble (left side of Fig. [2]) confirms this 
interpretation of the results. One clearly identifies a peak close to the center of the plane and the 
decrease towards the surface. No remnant of the building blocks, the merons, can be identified. 
The dynamical correlations change completely the picture (right side of Fig. |2]). The weight 
determined by the action in Eq. (|20l) favors configurations where destructive interference be- 
tween gauge fields of individual pseudoparticles prevent the formation of a large "mean field". 
In turn, the action density reflects clearly the structure of individual pseudoparticles. Screening 
of the gauge fields of the pseudoparticles is apparently also responsible for the restoration of 
translational symmetry. The comparison of the action density of merons and instantons in Fig.|3] 
demonstrates the similarity in the structure of the configurations generated by superimposing 
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Figure 2: Contour plot of the action density in the (x, y) plane for a single configuration of 
the stochastic (left) and of the dynamically correlated meron ensemble (right) with N M = 500. 
Left: s[A] = 110000, 15000 < s(x) < 327000. Right: s[A] = 1530, 180 < s(x) < 11300. 



and dynamically correlating meron and instanton fields, respectively. The larger amplitudes of 
the instanton fields give rise to correspondingly larger values of the action density. These results 
suggest decomposing the action density of the dynamically correlated ensembles into a single 
instanton (sj) or meron (sm) and a background (s B ) component. For instantons with their finite 
value of the action, the simplest Ansatz for the action density is [cf. Eqs. (1261271) 1 



The logarithmic infrared singularity of the meron action requires a slight modification in the 
Ansatz for the separation of single meron and mean-field contributions [fTTTl . The single meron 
contribution is obtained by integrating s(x) [Eq. ©] over a sphere of radius r, which for small 
meron size (p <C r) becomes 



g 2 {s) = g 2 ( Sl + s B ) = 8tt V + g 2 s B = 4.94 Nj + g 2 s B . 
For the dynamically correlated instanton ensembles of Table |3] we find 



(28) 




0.55. ..0.65. 



(29) 




(30) 



The matching requirement on r, 

s(r) = s B , 

yields the following expression for the action density, 




(31) 



For the dynamically correlated meron ensembles of Table [2] we find 



t4 = 0.65. ..0.75. 



(32) 
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Figure 3: The action density in the (x, y) plane for a single configuration of 200 merons (left) 
and 200 regular instantons (right) respectively. Averaged over the plane the action density is 
513 and 4700 respectively. 



3.2 Wilson Loops 

The following discussion focuses on the comparison of Wilson loops (flOl) in stochastic and 
dynamically correlated ensembles. A more thorough investigation of various issues concerning 
Wilson loops will be presented in later sections. Our standard choice for the path C in Eq. (flOl) 
is a rectangular path located in the (x,, Xj) planes with the center at the origin and with the ratio 
of the sides equal to 2. Thus, for a given configuration, we evaluate 12 different Wilson loops 
Wij and obtain our final results by taking the ensemble average and the average over the 12 
different orientations, 



is used to estimate the statistical uncertainty in the determination of (W). Fig. |4] shows the 
Wilson loops for the stochastic and correlated ensembles as function of the area. 

For sufficiently large loop size, we parameterize the logarithm of the Wilson loop as a sum of 
constant, perimeter (V) and area (^4) terms, 




(33) 



The variance 




(34) 



ln(W) = uj + tV - a A, 



(35) 



with 




for our standard choice of the path C. The parameters of the fit are 



u = 



-2.1, r 



2.3, a 



22.9 



(36) 
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Figure 4: Logarithm of a Wilson loop as a function of its area for three meron ensembles with 
Nm = 500: A stochastic ensemble (green x), an ensemble with dynamical color orientations 
(pink — ) and an ensemble with dynamical color orientations and dynamical meron positions 
(blue *). The fit curves are described in the text. 



for the stochastic ensemble, 

uj = -0.7, r = 1.1, a = 11.8 (37) 
for the ensemble with dynamical color orientations but random positions, and 

uj = -0.7, r = 1.1, a = 11.5 (38) 

for the ensemble with dynamical color orientations and positions. As in the case of the action 
density, the change from the stochastic to the dynamically correlated ensemble is significant. In 
comparison, only minor changes result when, in addition, the meron positions are also dynam- 
ical variables. For A > 0.1, the Wilson loops of the 3 ensembles are well reproduced by the 
parameterization d35l) . In particular, the results strongly suggest that all 3 ensembles give rise 
to an area law for sufficiently large size of the loops with a value of the string tension differing 
by a factor of 2 between stochastic and dynamically correlated ensembles. The large difference 
in the dimensionless ratios formed by action density and string tension, 

g 2 s cPs 

— — ~ 500 (stoch. ensemble) , — — ^11 (dyn. corr. ensemble), (39) 

<7 a 2 

reflects the essential difference in the dynamics of the ensembles. It is plausible that this strong 
decrease in the dimensionless ratio arises from suppression of fluctuations. One expects that 
reducing the fluctuations has a larger effect on the positive definite local action density than on 
the non-local Wilson loop. As will be seen below, the large value of the action density expressed 
in units defined by the string tension rules out the stochastic ensemble as a viable model of the 
Yang-Mills dynamics. 
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Small Wilson Loops in the Stochastic Ensemble 

In all the ensembles considered, the size dependence of small Wilson loops can be determined 
analytically. Since the configurations are built up by a finite number of pseudoparticles, fluctua- 
tions of the gauge fields cannot occur on arbitrarily small scales. In other words, for sufficiently 
small sizes, the gauge fields can be assumed to be spatially constant stochastic variables. For 
constant gauge fields, we use the identity 

e _ i<TC2 e ~i<rd e i,xC 2 e i<T Cl = 2i gin Ci gin (j^ e -kxC 2 ^ i x £^ a e i<rC 2 (4Q) 

and obtain with 

Ci = -giiAi 

W = 1-2 sin 2 X -gi x A x sin 2 X -gi 2 A 2 [1 - (A, A 2 ) 2 ] = 1 - t sin 2 X -gi x A x sin 2 X -gi 2 A 2 . (41) 

In the last step we have averaged over the color orientations. The space-time independent 
gauge fields Aj are generated by a superposition of a large number of the stochastic variables 
(positions and color orientations) specifying the individual pseudoparticles. We therefore expect 
the space-time components Aj to be Gaussian distributed, 

p(^)d 3 ^ = poA'dA^e-^^l 

We determine the width of the distribution by the expectation value of A 2 [cf. Eq. (l2~TI) l 

/ P (A)g^dAdn = 3a 2 = ^ ( } (42) 

and obtain 

Wo = 1 - I (l - (1 - iWo) e-« 2 ) (!_(!_ £ y Q ) ^Al^ , (43) 

with the values for the stochastic meron ensemble 

o"q = 120 for N P = N M = 500 , p = 0.16 . (44) 

The typical size of a gauge field component (with fixed color and space labels) is thus of the 
order of 

gA^l^/N^i^io. 

Depending on Np, to next order we still may consider ^-dependent fluctuations to be small as 
compared to Ci but to be of the order of 1 and therefore not suppressed. We write 

C f = Q + 7i « C f + Q 7 i , < |Ci|. 

To leading order, the directions of Ci and C« are identified but not their length. We find 

-tr e -19 ^ 2 e"^ 1 e if?<TC2 e i£,<TCl (45) 



COS71 cos 72 — ^Ci • C 2 j sin 7! sin 72 — 2 sinCi sin C\ ■ sin C 2 sin C 2 (l — (CiC 2 ^ 
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0.01 0.02 0.03 0.04 0.05 0.06 

Figure 5: Enlarged portion of Fig.|4]for small size Wilson loops. The curves show the result of 
the small loop expansion [Eq. (|46l) l (red, lower curve) for the stochastic ensemble and the result 
[Eqs. (I43I48I) ! (blue, upper curve) for the ensemble with dynamical meron positions and color 
orientations. 



By Taylor expansion of the gauge fields at the center of the cube, the fluctuations can be 
calculated. With 

we obtain for the deviations along a rectangular loop in the (1,2) plane 

Treating the deviations as independent stochastic variables with a Gaussian distribution 

p( t ) = far*'/** 

of variance 

2 _ I 2 \^ a2 

a 

we obtain 

(cos|) = (cos|)=e-^^ 2 / 4 . 

In this approximation, the expectation value of the Wilson loop [Eq. (1431) 1 acquires the following 
correction [cf. Eq. (|45T)1. 

W^Wo-l + e-^ 1 ^ 2 / 2 , (46) 
and we find for the stochastic meron ensemble 

o\ = 118 for N M = 500 , p = 0.16 . (47) 
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The agreement with the results of the numerical evaluation (cf. Fig. [5]) confirms the role of 
the constant gauge fields for small size Wilson loops. The above analysis of the small Wilson 
loops applies to some extent also to the dynamical ensembles where spatial variations of the 
gauge fields cannot happen either on arbitrarily small scales. Assuming independence of the 
fluctuation of the constant gauge fields, we still may use the expression (l43l) and, in the absence 
of a closed expression for the fluctuations of the gauge fields, determine the value of the width 

al = 10 (48) 

by a fit to our numerical results (cf. Fig. [5]). The reduction in a , i.e., in the size of the gauge 
field fluctuations, by one order of magnitude from the value (|44|) of the stochastic ensemble 
by the dynamics of the pseudoparticles is compatible with the reduction in the action density 
(dominated by the A 4 term), cf. Fig.Q] 



3.3 Thermodynamic Limit 

Unlike the stochastic ensemble, the dynamically correlated ensembles exhibit with the restora- 
tion of translational invariance a proper thermodynamic limit. In the case of stochastic ensem- 
bles, when doubling the linear dimensions of the system L with the density of pseudoparticles 
kept fixed, intensive quantities like the action density calculated in the stochastic ensemble 
do not remain invariant. The action density actually increases in this case by a factor 16 [cf. 
Eq. (|25l)l. Results for the dynamically correlated ensembles indicate significant improvement. 



S-2\ 



__. I i I i I i L 

'0 0,2 0,4 0,6 0,i 

A 



1 1,2 



1.4 



Figure 6: Logarithm of Wilson loop as function of the area for meron ensembles N M = 50 (L = 
2) (circle) and N M = 800 (L = 4) (star). 

The realization of translational symmetry (cf. Fig.Q]), or the scaling properties to be discussed in 
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the next section, would be hard to understand without a proper behavior in the thermodynamic 
limit. Here we present the results of a direct numerical examination of the thermodynamic limit 
by doubling the side length L and increasing the number of merons from 50 to 800. To achieve 



N M 


P 


<*> 


x 1/4 


C s (0) 


a-(o) 


50 


0.16 


273 


0.66 


2.50 • 10 5 


1.86 • 10 5 


800 (L = 4) 


0.16 


284 


0.66 


2.53 ■ 10 5 


1.88 • 10 5 



Table 1 : Test of thermodynamic limit for the vacuum expectation values of the action density s 
defined in © and of the topological susceptibility x (195b , and of the values of the correlation 
function of the action density (l8~8b and the topological charge density (189b at zero separation. 

high accuracy in this important test, we have generated ensembles which in comparison with 
our "standard" ensembles contain 4 times as many field configurations and where the action 
density in the update process has been evaluated at 3 times more meshpoints. The results in 
Fig. [6] and Table \T\ demonstrate the agreement in the Wilson loops and in various other observ- 
ables for these two ensembles. It is remarkable that the thermodynamic limit is established as a 
consequence of the dynamics of the pseudoparticles. 



4 Wilson Loops in Correlated Pseudoparticle Ensembles 

This section contains our main results concerning Wilson loops, an analysis of the structure of 
small Wilson loops and a discussion of the relation of our results to the general properties of 
Yang-Mills Wilson loops in the small and large loop limit as known from lattice gauge studies. 
The results have been obtained on the basis of about 150 000 - 250 000 configurations in each 
of the ensembles discussed in this section. 



4.1 Scaling Properties 

The Wilson loops evaluated in the meron and instanton ensembles with a wide range of pseu- 
doparticle numbers exhibit a universal behavior as demonstrated by Figs. [7] and [8j After rescal- 
ing the area A — > A A, the values of the Wilson loop lie on a universal scaling curve, where 
deviations from scaling are within the statistical uncertainties. To account for the rescaling of 
the area, we generalize the Ansatz (135b for the logarithm of the Wilson loop, 

\n(W) = lo + tVXV -aXA. (49) 

For merons, the values of the universal parameters are 

u = -0.95, r = 1.34, a = 12.7. (50) 

A central result of this work is that for both merons and instantons, for sufficiently large size, 
the logarithm of the Wilson loop decreases linearly with increasing area. This area law clearly 
demonstrates confinement in these pseudoparticle ensembles. The physical scale for each of 
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Figure 7: Logarithm of a Wilson loop as a function of its area for meron ensembles N M = 
800 (L = 4) (star), 100 (diamond), 200 (square), 500 (circle), 1000 (triangle). The values 
of the scale parameter A are given in Table [21 Also shown is the curve corresponding to the 
parameterization (|49l) with the values of the parameters given in (|50l) . 



the pseudoparticle ensembles is obtained by identifying the value of Act with the the value 4.4 
fm~ 2 , i.e., the unit of length (u.l.) for each of the ensembles is determined by A. For merons, 

lu.l. = 1.70 V\ fm. (51) 

The corresponding parameters for the instanton ensembles are 

u = -0.52, r = 1.18, a = 19.0, 1 u.l. = 2.08^^1 . (52) 

Fig.|9]displays the dependence of our results on the pseudoparticle size and coupling constant. 
Whereas changes in the Wilson loops induced by changes in the pseudoparticle size p can 
be accounted for to a large extent by a change in the physical scale, variations in the coupling 
constant by a factor of 4 modify the shape beyond changes in the scale, and a tendency is visible 
for a slightly smaller string tension in the ensembles with variable pseudoparticle size. We note 
that in the update process for the ensembles with variable size, the sizes of the pseudoparticles 
are kept fixed. On the basis of the identification of the value of the string tension extracted from 
the Wilson loops with the empirical value, we can compare our results for various observables 
with those obtained in other approaches, in particular with lattice gauge results. The values of 
the scaling parameter A together with other properties of meron and instanton ensembles are 
given in Tables [2] and [3j This first discussion of physical observables focuses on the action 
density and the topological susceptibility. Note that we refer in Tables |2] and [3] to results that 
will be explained in detail in section 5. From lattice SU(2) calculations for the action density the 
value (s)/a 2 = 4.5 lfi~4l has been deduced. QCD sum rule results Q~5l[T6]| range from (s)/a 2 = 
4.5 to 10. The necessity of subtracting divergent terms to define the continuum limit makes the 
uncertainty of these values of the action density, or equivalently of the gluon condensate, hard to 
assess. For instance, in the lattice calculations of Ref. [14], the divergent contributions are about 
a factor of 20-500 larger than the extracted value of (s) . As discussed above, the pseudoparticle 
action density also contains contributions which in the limit of vanishing pseudoparticle size 
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0,1 0,2 0,3 0,4 0,5 

A 

Figure 8: Logarithm of a Wilson loop as a function of its area for instanton ensembles Nj = 50 
(star), 100 (diamond), 200 (square), 500 (circle). The area has been rescaled with A given in 
Table [3j Also shown is the curve corresponding to the parameterization (1491 with the values of 
the parameters given in (l52i 



p — ► become singular and should therefore be subtracted. Thus a more relevant quantity to 
be compared with the lattice and sum rule results is the "mean-field value" of the action density 
[cf. Eqs. (1291321 1. We find 

— ~ 8 (10) for merons (instantons) . (53) 
<j 

This discussion also shows that the presence of undamped fluctuations in the field configurations 
of the stochastic ensemble [cf. Eq. (l39l) l is in severe conflict with the lattice and sum rule results. 
The topological susceptibility is directly related to the rf mass by the Veneziano-Witten formula 
lfi~7l[T8l and is well measured in lattice QCD. From the results in Tables [2] and [3] we conclude 

X 1/4 /^ 1/2 ~ 0.31 , f = l; 0.42 < x 1/4 /^ 1/2 < 0.48, £=2. 

Both the meron and the instanton ensembles yield values of the susceptibility of the correct 
order of magnitude. The value for the instanton ensembles agrees within its 7% error bars with 
the SU(2) lattice result EU x 1/4 /^ 1/2 ~ 0.48. 

4.2 Small Size Wilson Loops and the U(l) Limit 

Fields generated by superposition of a finite number of pseudoparticles are finite. Fluctuations 
do not occur with arbitrarily small wavelength. They are effectively cut off at length scales of the 
order of the pseudoparticle size or the average distance between neighboring pseudoparticles. 
As a consequence, the values of Wilson loops calculated in pseudoparticle ensembles vanish 
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Figure 9: Logarithm of a Wilson loop as a function of its area for meron ensembles Nm = 50. 
Left : Ensembles with different values of the meron size (0.01 < p < 0.16) and of the coupling 
constant (g 2 = 8, 32). The area has been rescaled with values of A in the range 0.31 < A < 0.44. 
Right: Ensembles of variable size uniformly distributed in the interval 0.03 < p < 0.17 (red) 
and of fixed size p = 0.10 (green). 
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[fm- 4 ] 


[fm- 1 ] 


[fm- 4 ] 


[MeV] 


1000 


2742 


1.25 


1.19 


3.67- 10 6 


2.49 • 10 6 


4.8 


0.30 


210 


118 


500 


1744 


1. 


1.01 


2.12 • 10 6 


1.45 • 10 6 


3.8 


0.27 


209 


120 


200 


897 


0.69 


0.9 


9.7- 10 5 


7.0 ■ 10 5 


3.2 


0.23 


226 


124 


100 


507 


0.52 


0.77 


4.8 ■ 10 5 


3.6 ■ 10 5 


2.8 


0.20 


225 


124 


50 


249 


0.37 


0.66 


2.23 ■ 10 5 


1.78 ■ 10 5 


2.7 


0.16 


218 


132 



Table 2: Properties of meron ensembles: Vacuum expectation values of the action density s 
defined in © and of the the topological susceptibility x (l95l ) (also in physical units), values 
of the scaling parameter A (l49l) , and values of the correlation functions of the action density 
(|88l ) and the topological charge density (|89l) at zero separation with the standard choice of the 
parameters (fT6l ) for ensembles of N M merons and with meron density hm = N M / L A . 



with vanishing loop size. This presence of an intrinsic cutoff prohibits a straightforward com- 
parison of Wilson loops of pseudoparticle ensembles with field theoretic treatments in which 
for small sizes, Wilson loops approach the limit of the Maxwell theory and yield Coulomb-like 
static quark-antiquark potentials that diverge with decreasing distance of the sources. Although 
both stochastic and dynamically correlated ensembles suffer from this deficiency, the stochastic 
ensemble with its much larger fluctuations appears to exhibit a Coulomb like behavior for areas 
A < 0.01 (cf. Figs. |4] and |5]) and only for values A < 0.002, does the Wilson loop show the 
quadratic dependence on the area. In the dynamically correlated ensembles we cannot identify 
such an intermediate Coulomb-like regime. For a meaningful interpretation we therefore have 
to account for the presence of this intrinsic ultraviolet cutoff and we will compare with Wilson 
loops in Maxwell theory in which the necessary ultraviolet regulator is kept finite. 
For small loop size we expand the exponential in the definition of the Wilson loop in SU(2) 
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5430 


1.0 


1.81 


2.7- 10 7 


1.6 • 10 7 


1.68 


0.33 


291 


162 


200 


2490 


0.66 


1.49 


9.7- 10 6 


6.6 • 10 6 


1.54 


0.27 


307 


164 


100 


1350 


0.48 


1.37 


5.1 • 10 6 


3.8 • 10 6 


1.45 


0.23 


314 


180 


50 


651 


0.32 


1.15 


2.3 • 10 6 


1.9 • 10 6 


1.64 


0.19 


340 


190 



Table 3: Properties of instanton ensembles, as defined in Tabled 
Yang-Mills theory 

W = itr jpexpi^ ^A^xjdaA « 1 - ^ j> A^dx" j> A^dy" . (54) 

Here, A denotes one of the color components of the gauge field. For small loop size we compute 
the vacuum expectation value of the Wilson loop by identifying A with an abelian gauge field 
and find 2 

(W)ttl-^Lj J cu (x)K(x-y).r c (y)d 4 xd 4 y, (55) 

where 

JN X ) = f 8\x-x c (s)) — ds 
J ds 

is the current of a charge transported along the curve C. The Euclidean gauge field propagator 
K in Lorentz gauge is given by the inverse of the 4-dimensional Laplacian. The integral (1551) 
requires regularization. We choose a heat kernel regularization and obtain for the regularized 
propagator in Maxwell theory 



K( x ) = / e 



a 



(56) 



(2tt) 4 k 2 An 2 x 2 
We consider rectangular loops with side lengths R and T and denote with 

uj{R,T)= [ dx [ dy K(./(x -y ) 2 + R 2 )) (57) 
Jo Jo 

the contribution to W of the time-components of the two currents on opposite sides of the 
rectangle in the expression (|55T ). The Wilson loop is obtained from the 4 contributions of time 
and space components on the same and on opposite sides of the rectangle, 

3o 2 

(W) w 1 - —[u(0,T)+u(0,R)-u(R,T)-u(T,R)]. (58) 



The integral (1571 ) can be evaluated in closed form with the result 



u(R,T)= 1 



2tt 2 



, a + Va 2 + R 2 + T 2 T ^ ^ VR 2 + T 2 + a 2 -a 

In , arctan T R 



a + Va 2 + R 2 R R 2 ^R 2 + T 2 + a 2 + aT 2 - 

(59) 

Performing the limit a — » 0, we obtain for R < T 
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Figure 10: Logarithm of rectangular Wilson loops of fixed area (A = 0.01, 0.02, 0.04, 0.06) as 
a function of the perimeter for an ensemble with 500 instantons of variable size p < 0.1. The 
curves are obtained from Eqs. (I58ll59l) with the value (l62l) of the regulator. 



(W) « 1 + 



3g 2 tttT 



— + 1- 



R + T 




(60) 



16tt 2 12R 



a 



which contains both the Coulomb potential and the self-energy of the static charges. On the 
other hand, in the limit of small loops and for fixed a, (W) — 1 tends to zero with the square of 
the area 



as is the case for Wilson loops in pseudop article ensembles. The small Wilson loops in the 
stochastic ensemble in Fig. \5\ display these two different limits. The regime (|6"TT) of vanishing 
loop size is restricted to values of the area of the order of 10~ 3 and is followed by the Coulomb- 
like behavior associated with the limit (l60l) of vanishing regulator size, which extends up to 
values of the area of the order 10~ 2 . With the reduction in the fluctuations, the dynamically 
correlated ensembles do not explicitly show the Coulomb-like behavior. Nevertheless on the 
basis of this analysis we can understand the small loop limit of the dynamically correlated 
ensembles and make the connection to the perturbative U(l) limit. 

For the comparison of the pseudoparticle Wilson loops with Wilson loops of the Maxwell theory 
(1581591) at finite a, we have determined the value of the regulator by a fit of expression (|58l) to 
the Wilson loops in the pseudoparticle ensembles. The results are shown in Fig.QTJJ The value 
of the regulator is 



The heat kernel regulated expression for the Wilson loop (1581) reproduces the trend of the nu- 
merical results. Quantitative agreement is obtained for large values of the perimeter, i.e., for 
R < T. 



4.3 Large Size Wilson Loops and the String Limit 

Results from lattice gauge theories (cf. e.g. Il20l ) support the conjecture that large Wilson loops 
can be described by a an effective string theory. In the string picture, one assumes the Wilson 




(61) 



a = 0.15. 



(62) 
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loop is determined by the string partition function ETTl 

W{R, T) = I d[x ± }e-^ RT - s ^\ S qf = [ R f d& d& d a xl d a xl . (63) 
J Jo Jo 

Apart from the classical string energy aR, this expression accounts for the small quantum fluc- 
tuations around the classical string. The integral over the fluctuating transverse string coordi- 
nates £^(£1, £2) is given by the determinant of the two-dimensional Laplacian, 

e- = /dMe^» = aif ^L_. (64) 

For computation of det(— A), one imposes Dirichlet boundary conditions along the boundary 
of the Wilson loop. In terms of the resulting spectrum, the free energy F for a rectangular loop 
of side lengths R, T reads 

m,n=l x ' 

In C function regularization [|22l|. the following finite expression 

det(-A) = — =r](i^] (66) 



.f2R \ R 

in terms of the Dedekind 77 function 



2mnz\ 

L — e 

n=l 



is obtained. The prefactor in the i] function yields the well known Liischer term ll20ll . ttT/12R, 
which dominates the logarithm of the Wilson loop in the limit of large T/R. 

For the purpose of analyzing the Wilson loops in pseudoparticle ensembles, we need a regu- 
larization scheme that makes explicit the role of the cutoff. In heat kernel regularization, the 
expression for the free energy reads 

m ^ , / V 2 m 2 7T 2 n 2 \ > / .y 1 ^ 

TF(\) = ln 1-fT + -Je ) e- A V^- + ^ . (67) 

m,n=l ^ ' 

The leading term in A -1 can be evaluated in closed form 

TF(X) « ^ [dpdq \n(p 2 + q 2 ) e~ A v / ^+^ 

IT 2 J 

ATR( d\ [°° _ A 4TR 1 , . 



In the limit A — > 0, a quadratically divergent correction to the string tension is obtained. As 
follows from dimensional arguments, the subleading linearly divergent (up to logarithms) terms 
are multiplying R and T. 

For comparison with the large Wilson loops evaluated in the pseudoparticle ensembles, we 
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Figure 1 1 : Logarithm of rectangular Wilson loops of fixed area A as a function of the perimeter 
for an ensemble with 500 instantons of variable size p < 0.1. The 7 values of the area vary 
between 0.08 and 0.40. The solid and dashed curves are obtained from Eqs. d70l) and (|68l) with 
the values of the parameters given in Eqs. (1711) and (l69l) respectively. 



keep A finite and introduce counterterms associated with the divergent contributions. We fit our 
numerical results for the Wilson loop with the Ansatz 

In W = -TF(X) + a + b(T + R) + a RT , (68) 

with the counterterms a, b and a depending on A. We have applied this procedure to Wilson 
loops of ensembles with 500 instantons and merons respectively with variable size (p < 0.1). 
The fit to the pseudoparticle results yields the following values of the parameters, 

A = 0.25, a = -0.7, b = 1.8, cr = -14 . (69) 

As Fig. [H] shows, the string model defined by Eqs. (1681691) leads essentially to a linear de- 
pendence of the logarithm of the Wilson loops as a function of the perimeter with a very mild 
dependence of the slope on the value of the area. For sufficiently large aspect ratios of the loops, 
the string model reproduces the numerical results for the pseudoparticle ensembles. However 
the model misses the change in slope close to the threshold, i.e., in the limiting case of a square. 
For fixed area A, the slope of the Wilson loop as a function of the perimeter V becomes in 
general infinite if the Wilson loop depends on |T — R\. We find 

i.e., a divergent slope at threshold is obtained for [i < 2. The presence of such a term is indicated 
by the rapid change of the numerical results with V at the threshold V = A\f~A. To account for 
this structure, we modify minimally our Ansatz (l68~l) . 



In W = -TF{\) + a + bVT 2 - R 2 + a RT, (70) 



25 



Figure 12: Left: Logarithm of rectangular Wilson loops with aspect ratio 1:2 as a function of 
the area for an ensemble with 500 instantons of variable size p < 0.1. The solid curve has been 
obtained from a fit of In W in Eq. (|7Q|) to the data in Fig. [TT] with the resulting values given in 
(1711) while the dashed curve shows our standard fit (|74l) to the data points with aspect ratio 1:2 
in the present figure with the extracted values (1751) . Right: Asymptotic extraction of the string 
constant from Wilson loops as a function of the area. Notice the much larger scale of the area. 
Crosses represent the function (|7Q|) with the values of the parameters (fTTT) . The straight line is 
obtained from fitting the asymptotics of (1701) (crosses). The extracted value of the string tension 
is given in Eq. (|76|) . 



such that neither the asymptotic behavior at large area nor, for fixed area, at large perimeter is 
changed. In turn, the changes in fit parameters are very small [cf. Eq. (|7T| )1 



A = 0.25, a = -0.7, b = 1.8, a = -13 . (71) 

As Fig. [TT1 demonstrates, the Ansatz (1701) catches the essential properties of the pseudoparticle 
Wilson loops of sufficiently large size. Deviations at large values of the perimeter must be 
expected. The string picture breaks down if the smaller length R of the rectangle is of the order 
of the diameter of the flux tube, 

R « 2r , 

and therefore the perimeter has to satisfy 

V < V max = . (72) 

r 

With r ~ 0.15 fm (corresponding to 0.07 in our units) we obtain for A = 0.15 (second curve 
from the top in Fig. fTTj) V m£LX = 2.4, a value which, as can be seen from Fig.QT] is of the correct 
order of magnitude. 

The origin of the observed threshold behavior is not obvious. It could be a peculiarity of the 
effective regularization in pseudoparticle ensembles. Within the string model description, to 
lowest order, the string partition function [Eq. ((671) 1 does not depend on the variable \T — R\. 
One can show that the derivative of the ( function (l66l) with respect to the perimeter is finite at 
threshold. Infinite slopes of the Wilson loops however appear in higher order corrections which 
have been calculated in £ function regularization ll22l . It would be interesting to extend our 
analysis and include such corrections within heat kernel regularization. 
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Figure 13: Logarithm of rectangular Wilson loops of fixed area as a function of the perimeter for 
an ensemble with 500 merons of variable size p < 0.1. Left: A = 0.005, 0.01, 0.02, 0.04, 0.06. 
The curves are obtained from Eqs. (I58~ll59~l) with the value (|77l) of the regulator. Right: The 10 
values of the area vary between 0.08 and 0.72. The dashed and solid curves are obtained from 
Eqs. (l68l and (|70l) with the values of the parameters given in Eqs. (1781) and (|79|) . respectively. 



Having determined the counterterms (fTl!) . we extract the value of the string constant a from the 
asymptotic behavior of In W [Eq. (|7Q|) 1 

aRT = lim -T F(\) + a RT . (73) 

RT->oo 

Fig. d2l illustrates the procedure. On the left hand side our standard procedure for determining 
the string constant by fitting the Wilson loop at constant aspect ratio (1:2) is shown. The fit with 
the standard Ansatz 

In W = lu + rV - a A (74) 

and extracted values 

u = -0.1, r = 0.58, a = 20.5 (75) 

yields an excellent fit to our numerical results. The extracted value of the coefficient r is about 
a factor two smaller than the slopes of the dashed curves in Fig. [TH and does not adequately 
describe the numerical results. 

The Ansatz based on the string model ([70]) . which also accounts for the dependence on the 
perimeter for fixed area, describes the pseudoparticle Wilson loops in a limited region only. 
Clearly the string model is not suited to describe the perturbative region of small loops. On the 
right hand side we show the extraction of the string constant from the asymptotics. Within the 
uncertainty of the procedures, the value 

a = 21 (76) 

agrees with the value in (1751) . 

Very similar results are obtained for the Wilson loops in meron ensembles. In Fig. \T3\ the 
results of the analysis are presented for small and large Wilson loops respectively. The value of 
the heat kernel regulator for the small Wilson loops [cf. Eqs. (1581591 )1 is 

a = 0.20 . (77) 
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The parameters for the string model Ansatz (1681 have been determined to be 



A = 0.21, a = -1.2, b = 0.9, a = 0.8 



(78) 



and for the modified Ansatz (1701 



A 



0.21, a 



-1.1, 6 = 0.9, a = 1.3. 



(79) 



The string tension extracted from the asymptotics of Eq. (1701 is 



a = 12.5, 



(80) 



in agreement with the value deduced from the parameterization (1741) 



to = —0.5, t 



0.94, a 



12.8. 



(81) 



The consistency of the results of these two different methods is an important indicator of the 
robustness of our determination of the string tension. 

4.4 Wilson Loop Distributions and Higher Representation Wilson Loops 

Wilson loop distributions and Wilson loops in higher representations further characterize the 
confining Yang-Mills dynamics. In particular, by lattice calculations, an intermediate regime 
of Casimir scaling has been established for both SU(2) and SU(3) Yang-Mills theories and 
deviations from Casimir scaling by string breaking have been identified, cf. Il23l l24l l25l l26ll . 
More recently, the Wilson loop distributions have been introduced lfi~3l as quantities that contain 
the full information on the Wilson loops in higher representations. The study of Wilson loop 
distributions revealed an unexpectedly simple property of the Yang-Mills dynamics. These 
distributions can be described as a result of a diffusion process. Here we will present the Wilson 
loop distribution for pseudoparticle ensembles and discuss the related expectation values of 
Wilson loops in higher representations. We will compare our results with the diffusion model 
for the Wilson loop distribution. In this model the Wilson loop distribution p(cos#, t) is a 
solution of the diffusion equation on the group manifold S 3 , 



The diffusion model makes no assumption about the connection between the "time" t and the 
Wilson loop size. For given size of the loop, we determine the value of the time t by fitting the 
expectation value of the Wilson loop. As demonstrated in the left part of Fig. El the Wilson loop 
distributions are well described by the diffusion model. In particular, the change in shape from 
the strongly peaked distribution for small loops to the Haar measure for large loops is correctly 
reproduced. These results are very close to results obtained in lattice gauge calculations lfi~3l 
and, together with our other results, strongly support the treatment of the Wilson loop dynamics 
in the pseudoparticle approach. 




(82) 



and is given by 




(83) 



71=1 
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Figure 14: Left: Wilson loop distribution for an ensemble of configurations containing 500 
merons. The sizes of the loops are 0.48 x 0.24,0.72 x 0.36, 1. x 0.5 with average Wilson 
loop values 0.66, 0.26, 0.04. These results are fitted with the distribution p(cost9, t) (f8~3l) with 
parameters t = 0.13, 0.44, 1.05 respectively. Right: Variance of the Wilson loop Aw (l8~4l) as a 
function of the average (w) for ensembles of merons with Nm = 50, 100, 200, 500, 1000, 3200 
together with the results of the diffusion model Eq. (l8~4l) . 



Given the distribution of Wilson loops, the variance of Wilson loops is easily computed, 

Aw = ^/(w 2 ) - (w) 2 = \\J^~ 4 W 2 + 3(w) 8 / 3 , (84) 
as well as the expectation value 

(Wj) = e~^ j+1)t (85) 
of Wilson loops in the (2j + l)-dimensional representation, 

Wm = — tr exp{2ztf (~ j ■ • )} = f^ 3 ± ^ . (86) 

The right hand side of Fig. [T4l demonstrates the universality of the variance of the Wilson loop 
for meron ensembles with meron numbers Nm varying by up to a factor 60. Agreement with 
the diffusion model within 10% or less is found. The numerical results for Wilson loop ex- 
pectation values in different representations are shown in Fig. [131 together with fits using the 
parameterization (|49l . With increasing dimensionality of the representation, the negative slope 
of the logarithm of the Wilson loop expectation values increases by about a factor 5 (cf. Table 

E3). 



J 


UJ 


T 


a 


1/2 


-0.72 


1.10 


11.5 


1 


-0.99 


1.89 


24.5 


3/2 


-0.85 


2.24 


37.5 


2 


-0.96 


2.95 


54.5 



Table 4: Parameters for the fit (|49l) of Wilson loops in the representation j for a meron ensemble 
with N M = 500 (A = 1). 
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Figure 15: Logarithm of a Wilson loop as a function of its area for meron ensemble with 
iVjif = 500 in j = 1/2,1,3/2, and 2 representations. Also shown are curves corresponding to 
the parameterizations (|49|) with the values of the parameters given in Table HI 



According to the diffusion model Eq. (1851) . Wilson loops in different representations exhibit 
Casimir scaling, i.e., the values of the string tension, or more generally the interaction energy 
of static charges in higher representations, are proportional to j(j + 1). We have 



R i = ,_/„m = ttz ■ (87) 



HW j+ i) = (j + l/2)(j + 3/2) 
\n(Wj) + 

Casimir scaling of the Wilson loops in two meron ensembles is examined in Fig. [16] On the 
average, we find that Casimir scaling is satisfied within 8%. We observe a decrease in Rj 
by about this amount. It appears that Casimir scaling is exact for small Wilson loops. Also 
these results are compatible with the findings in lattice gauge calculations where the validity 
of Casimir scaling has been demonstrated for SU(2) ll2"TI and SU(3) [|23l Yang-Mills theory. 
Indications for deviations from Casimir scaling for sufficiently large loops have been obtained 
in [|24l and [|25l and string breaking has been observed in [|26l . Due to poor statistics in the 
computation of large loops, we have not been able to establish or rule out string breaking in the 
pseudoparticle ensembles. 

5 Correlation Functions 



Correlation functions are useful tools to study susceptibilities, response functions, and excita- 
tions of many particle systems and quantum field theories. In the case of Euclidean time, they 
also enable the measurement of the energies of low-lying excited states. In this section, we 
will use appropriately chosen operators to study the topological susceptibility and low glueball 
masses of our effective theory, and compare them with corresponding lattice results. We be- 
gin with a discussion of methodology, describe the operators used to probe the system, discuss 
several related theoretical issues, present results for the + , 1 + , 1 _ , 2 + , and 2~ glueballs, and 
discuss a correlation function that reflects confinement. 
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Figure 16: Ratios of logarithms of Wilson loops of merons in consecutive representations Rj 
[Eq. (|S7H for j = 1/2, 1, 3/2, N M = 500 (left) and N M = 200 (right). 

5.1 Correlation Functions of the Action- and Topological Charge Density 

As an introduction, it is useful to begin by considering two correlation functions that we can 
analyze to a large extent analytically, the correlation functions of the action density [Eq. ©], 

Cs (x) = ([s(x)-(s)][s(0)-(s)]), (88) 

and of the topological charge density [Eq. ©], 

C s (x) = (s(x)s(O)). (89) 

The action density s(x) and topological charge density s(x) are the simplest examples of gauge 
invariant quantities. They are scalar or pseudoscalar quantities and bilinear in the field-strength. 
The explicit expressions of these quantities for a single meron or instanton are given in Eqs. © 
and ©. The normalization of the correlation functions yields the variance of the action and 
topological density respectively. Many results from this section have been tabulated in Tables |2] 
and|3]for previous reference and discussion, and one can see from these tables that fluctuations 
of the action are of the same order as the average value, 

The large fluctuations are due to the large differences in the values of the action density inside 
the pseudoparticles and in the background. In the center, 



s (0) = s (0) = 




which, in the ensemble of fields with 500 merons for example, is an order of magnitude larger 
than the average action density. The square of the action and of the topological density is 
concentrated in the core of the pseudoparticles. Therefore the strength of the fluctuations and 
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more generally the correlation function can be estimated by the incoherent sum over the sin- 
gle pseudoparticles. In such an incoherent single pseudoparticle approximation, we define the 
correlation function of, for instance, the topological density by 



N P 



1 / NT 

cr°\x) = V J2 ~ s ^y - *) s °(y + x - *) = -f \ d4y So(y) 5o(y + x) > (90) 
i=i j •* 

and similarly for any other observable. In this approximation, the integration over the color 
orientation gives a multiplicative factor which is canceled by the normalization of the partition 
function. Using the explicit expressions © and ([8]), we obtain the following estimates for the 
fluctuations in meron ensembles (N = Nm) 

^incoh/nx _ 155l7r2 N M Wncoh/ n X _ 17 ' ^ N M fQU 

Cs (0) -^80~7v" ' (0) "^5~^' (91) 



Cr° h (0) = CT°\0) = — £ = 2. (92) 



and for instantons 

384tt 2 Nj 
7 p A V 

Comparison with the numerical results in Tables [2] and [3] shows that the general trends in the 
values of the correlation functions for vanishing separation are correctly reproduced by this 
incoherent superposition of pseudoparticles. The similar values of C s (0) and C s (0), the dif- 
ference by an order of magnitude for meron and instanton ensembles, the dependence on the 
pseudoparticle size and number density agree typically within 20-30 % with the numerical re- 
sults. Although the pseudoparticle ensembles are far from being random, the fluctuations are 
dominated by the variations of the densities on the small scale of the meron size where the 
longer range correlations between pseudoparticles have little effect. 

In the following, we will compare our numerical results for the correlation function with the 
single pseudoparticle predictions [cf. Eq. (|90l )l which we have evaluated numerically for various 
observables. Besides the values for zero separation [Eqs. (|9Tl[92l l, the asymptotic behavior can 
also be determined analytically. We find for meron ensembles 

+ ^mcoh/.N _^ 971-2 N m \nt/p nincoh{ , oa 2 N m P 2 c . rcm 
t^oo. C s [t)^—— C s (t)^ 247r— — , £=1, (93) 

and for instanton ensembles 

AT 4 

t -> oo : Cl ncoh (t) = Cr oh {t) -> 24tt 2 -^^ , i = 2. (94) 

The integral over the topological charge density correlation function yields the topological sus- 
ceptibility 

X= (A)' / d 4 xC s (x), (95) 



S7T 



an important quantity of the strong interaction. In Tables [2] and [3] the values of % 1/4 for the 
various ensembles are given. Once more, these values essentially arise from an incoherent sum 
of single pseudoparticle contributions. For the calculation of x we need the curvature /3, 
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Figure 17: The correlation function \C s (x)\ for Nm = 50, p = 0.01, 0.02, 0.04, 0.08, 0.16 in 
comparison with the single-meron correlation functions. Note that in every case the correlation 
function has a node beyond which it is negative. 

and by approximating the single pseudoparticle correlation function by a Gaussian of the same 
curvature we obtain 

*' /4 -& (0)1/4 

resulting in 

/? incoh = 184J_ incoh= 4 

177 p" 1 3 p^ 

and the following expressions for meron and instanton ensembles 

/ . ,\l/4 / \l/4 / . ,\l/4 / n1/4 

(x mc ° h ) = 0.52(N M /V) , ^ = 1; (x mcoh ) = 0.83[N I /V) , £ = 2. 

(97) 

The iV 1 / 4 dependence is easily understood on rather general grounds. If we assume the topo- 
logical charge in volume V to be given by the difference in the number of pseudoparticles and 
anti-pseudoparticles N — N, the fluctuation in the charge is expected to be ~ \/~Nm an d hence 
X ~ N. The estimate d97l) deviates from the numerical results by about 10 %. For small 
pseudoparticle sizes, the approximation of the distribution by a Gaussian becomes invalid and 
deviations of the order of 20-30 % occur. 

The case of the topological charge density correlation function is particularly interesting. Since 
s ~ E • B, it acquires an i in Euclidean space and the correlation function is negative for 
nonzero x, Cs(x) ~ — J2n I (0|^l n ) \ 2 e~ EnX . The positive topological susceptibility arises from 
a positive contact term at at the origin [|27l . As Fig. [FT] shows, the pseudoparticle approximation 
to the path-integral exhibits these properties. The contact term is approximated by the single 
pseudoparticle distribution [Eq. (|90ll, and after the change of sign, the correlation function 
exhibits a new length scale. For p = .02, the results follow the single meron correlation function 
until s(x) has decreased to a value of about 5 ■ 10 -5 . Only at this level does a new length scale 
become important. We see here in detail how with decreasing cutoff p, the strength remains in 
the short range peak to produce a topological susceptibility essentially independent of the value 
of the cutoff. 
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Figure 18: Emergence of the hadronic scale in the correlation function C s (x) in ensembles 
with 50 instantons of size p = 0.05 (left) and 800 merons (L = 4) of size p = 0.2 (right), in 
comparison with the respective single pseudoparticle correlations functions. 

In section 4, we have already compared the measured topological susceptibilities for instanton 
and meron ensembles with lattice results. The instanton result, 0.42 < x 1 ^ 4 /^" 1 ^ 2 < 0.48, agrees 
within roughly 10% with the SU(2) lattice result lfi~9~ll , whereas the meron result x 1 ^/ (j1 ^ 2 ~ 
0.31 is significantly lower. From the developments in this section, we now see the origin of this 
difference. The estimate (1971) of the ratio of 1.6 of the topological susceptibilities (x 1 ^ 4 ) of in- 
stantons and merons arises from the fact that the topological charge of instantons is twice that of 
merons, and this ratio is only slightly reduced when expressed in physical units. Fig. [TBI shows 
results for the action density correlation function that illustrate the emergence of a hadronic 
scale, that is, a scale generated by dynamics that is significantly different from the scale set 
by the size of the pseudoparticles. If an effective theory has the correct degrees of freedom, 
we should expect it to be able to generate scales larger or smaller than that of the underly- 
ing pseudoparticles. The left figure demonstrates the superposition of instantons of small size 
(p = 0.05 fm) to generate a much larger correlation length, a familiar phenomenon arising in 
many-body systems. The right figure shows an example of the opposite case, where superim- 
posing merons of size 0.2 fm actually generates a hadronic scale that is smaller than the meron 
size. To emphasize the very slow decay [cf. Eq. ( |93l) l of the single meron correlation function, 
its profile is also shown in the figure. The emergence of the hadronic scale is a rather sub- 
tle effect generated by destructive interference among the pseudoparticles. To demonstrate the 
dominance of this hadronic scale over a large distance despite the weak decay of the building 
blocks, we had to use a meron ensemble in a volume 16 times larger than usual and containing 
16 times more pseudoparticles. 

The successful generation of the negative correlation function for the topological charge density 
and of hadronic scales both larger and smaller than the pseudoparticle scales for the action den- 
sity are important indicators of the potential of the pseudoparticle effective theories to describe 
low energy dynamics. Based on this success, we now use correlation functions to explore the 
low lying glueball spectrum. 
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5.2 Correlation Functions and Glueball Masses 



To the extent to which the effective degrees of freedom succeed in representing the full path 
integral, the integration yields a transfer matrix and one may write a Euclidean correlation 
function as a sum over eigenstates, 

(O(x)O(0)) ~ J2( Q \°\ n ) e ~ EnX ( n \°\ Q ) ■ ( 9g ) 

n 

As in lattice gauge theories, glueball masses may be extracted from the exponential decay at 
sufficiently large distances. In the case of the action density, we thus expect for large separation 
x that C s (x) ~ | (ra |s|0) | 2 e~ moa: , where m is the mass of the lowest ++ glueball. Similarly, 
the asymptotic decay of the pseudoscalar correlation function is expected to be determined by 
the lowest CT + glueball mass. 

5.2.1 Methodology 

To facilitate the extraction of single exponential terms in Eq. (|981) , it is advantageous to project 
onto states with zero momentum so that the continuum of states of non-zero momentum as- 
sociated with each intrinsic excited state is eliminated, and the gap between the desired state 
and the next higher energy state is as large as possible. In addition, although in a rigorously 
translationally invariant system it would be sufficient to project either the source operator at the 
time origin or the sink operator at t, since the current calculation is only approximately transla- 
tionally invariant, it is beneficial to project both the source and the sink. Hence, we will use the 
doubly momentum projected correlation functions 

C (t)= [ d 3 x [ d 3 y<O|0(*,x)0(O,y)|O). (99) 
Jv Jv 

The time direction in 4-dimensional space is defined here by the coordinate over which one 
does not integrate. The integration over the 3-dimensional cube V is limited to i « 0.5L...0.7L 
in order to avoid surface effects. 

In order to accurately determine glueball masses, it is necessary to identify a substantial region 
of pure exponential decay between short distance artifacts associated with the intrinsic size p of 
our effective degrees of freedom and large distance artifacts associated with the finite volume 
edge effects and the fact that beyond some point the fluctuations arising from finite statistics 
are comparable with the magnitude of the decaying exponential. In general, the region between 
the short distance and large distance artifacts is quite limited, and the identification of the ex- 
ponential regime is ultimately subjective. Unfortunately, this precludes quantifying statistical 
and systematic errors, and in the following discussion we will simply explain our rationale for 
selecting the fiducial regions we use to extract the slopes defining the exponential decay of each 
state, and give representative values, or ranges of values for each slope. Fig.[T8]is a good exam- 
ple of a case in which there is clear separation between the short range artifacts associated with 
the intrinsic size p and a substantial region of nearly exponential decay. Similarly, in Fig. \T7\ 
as previously noted, one observes the expected negativity in C s (x) beyond 2p, indicating that 
one is beyond the range of short distance artifacts. However, in this case, it is also clear that the 
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fluctuations become so large at larger distances that there is no correspondingly clean region 
of exponential decay, and we are not able to extract a mass from the pseudoscalar correlation 
functions in the present work. 

Given the difficulties in determining masses in many cases, we have also required that we 
obtain consistent results (< ±10%) for ensembles with different meron sizes or ensembles 
with different number of pseudoparticles. When possible, we have also used different operators 
to create states of the same quantum numbers, and sought consistency between these results. 
Finally, for spin 1 and 2 excitations we required deviations from rotational invariance to be 
small in the fiducial region. 

5.2.2 Operators 

To calculate glueball masses, we use correlation functions based on two kinds of operators: 
local operators constructed from products of two or three field-strength operators, and non- 
local Wilson loop operators. They are closely related to gluonic operators used in recent lattice 
calculations of the glueball spectrum and matrix elements (see ll28l for the relationship to the 
lattice symmetry group and further refinements). We first introduce the local operators and 
discuss their behavior under interchange of Euclidean space and time axes, and then describe 
the Wilson loop operators. Properties of the operators and calculations based on them are 
summarized in Table [5J 

Local operators with two or three field-strengths 

We have considered the complete set of gauge invariant operators consisting of products of 2 or 
3 field-strength operators, which we represent in terms of magnetic and electric fields in order 
to make the transformation properties under 3-rotations explicit. We use the following set of 
observables consisting of two field operators 

s = ^(E 2 + B 2 ) , T+ = E<g>E + B<g>B, T_ = E<g>E — B<g>B , T a = E®B + B®E, (100) 
with the traceless tensors 

(X ® Y) y = X?Y„ a - lx. a Y a 5 l3 , (101) 

6 

and the following observables consisting of products of 3 field- strengths 

V p = e abc (E a ■ B b ) B c , V a = e abc (E a • B b ) E c , 
T 3 = e*[(B a xB i )8E c + E c ®(B"xB fc )]. (102) 

To improve the statistics we have summed the correlation functions for vectors and tensors over 
the spatial components, e.g. 

3 

((X®Y)(a;)(X<8>Y)(y))= £ ((X ® Y) y (a;)(X ® Y) -(y)) 

j>i=i 

and have used the equality of the components as a test for rotational invariance of the ensembles. 
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Our subsequent calculations all refer to instanton ensembles, for which the fiducial interval 
for extracting masses is larger. In cases where masses could be extracted from both kinds of 
ensembles, the results agreed within the uncertainty of the procedure. We observe that the 
instanton fields [Eqs. ([5l[6])] satisfy 

B t = = -4 7-2^)2^ • 

As a consequence, gauge invariant combinations consisting of products of two or three field- 
strength operators have to be either scalar or pseudoscalar. For merons this is only true after 
momentum projection. If the momentum projection is not complete, the presence of single 
meron contributions will complicate the analysis of the correlation functions. 

In addition to the action density and the three tensors [cf. Eq. (11001) 1. we can form a sec- 
ond scalar, B 2 — E 2 , the Hamiltonian density transformed to Euclidean space, and the Poynting 
vector, E x B. Although if integrated over all space, the vacuum is an eigenstate of the resulting 
operators, our ensembles of pseudoparticles do not yield exact eigenstates of the Hamiltonian 
and projection to zero momentum is not exact. Hence, we obtain non-vanishing correlation 
functions. However, we find that the fluctuations in the Hamiltonian, i.e., the value of the cor- 
relation function for vanishing \t\, are one to two orders of magnitude smaller than those of the 
action, depending on the ensemble. In the case of the integrated Poynting vector, the momen- 
tum operator, we find very strong dependence of the correlation functions on the parameters of 
the ensembles, which highlights the unphysical origin of the correlations. For these reasons we 
do not consider these correlations further. Finally we remark that the "singular" behavior (cf. 
Fig.fTTI) of the pseudoscalar (E • B) correlation function C$ prevents a meaningful extraction of 
a mass after momentum projection. 

Behavior under interchange of Euclidean space and time axes 

In order to understand some of the results obtained subsequently, it is useful to note that by 
relativistic covariance, electric and magnetic fields can be transformed into each other. For van- 
ishing separation, the distinction between electric and magnetic field depends on the choice of 
the coordinate system. We illustrate the equivalence by considering the following contribution 
to the correlation function of T a 

(%(ar 4 )) = (A£j(0,x 4 )A£,(0)) = (F 4i F 4j (0,x 4 )F 4i F 4i (0)) . 

Here color indices are suppressed and no summation over i and j is carried out. Under a 
transformation which rotates the x 4 and X{ coordinates into each other 

. CC{ ^ %i ~~ ^4 5 ^4 ^ ^4 ~~ (103) 

we find, assuming 4-dimensional rotational invariance in the ensembles 

R : (t l3 {x 4 )) = (F 4i F l3 {x') F 4i i^(0)> , x' k = 5 ki x 4 , i^j. (104) 
For vanishing separation (x 4 = 0), this result, if expressed by electric and magnetic fields 

3 

(5^,(0) BtE^O)) = e ijk (BiB k (0) B % B k {Q)) , i^j. (105) 

k=l 
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together with the definition (1 1001) of the relevant tensors is seen to yield the identity 

C r+ (0) = C Ta (0) (106) 

for the off-diagonal elements. The diagonal elements of the traceless tensor can be connected 
to the off diagonal ones by rotations not involving the 4-coordinates. We note that the rotation 
(11031) has no counterpart in Minkowski space. This method is not directly applicable for non- 
vanishing separations since electric and magnetic fields are defined with respect to the direction 
of separation chosen as time axis 

Here z denotes the 4- vector defining the separation and z the corresponding unit vector. Only 
with this definition electric and magnetic field do not mix under 3-dimensional rotations in the 
space transverse to the direction of separation, i.e., under rotations satisfying 

R)jlvZ v Z^ . 

In this way the resulting operators project on states with definite angular momentum and par- 
ity. This redefinition invalidates the above derivation for finite separations. More complicated 
relations between the two correlation functions can be derived which however defy a direct 
numerical application. As we discuss below, it is striking that our numerical results for the 
two momentum projected correlation functions Ct + (x) and Cx a (x) exhibit essentially identical 
normalizations and a strong similarity in shape. Significant deviations occur only at the level 
of about 0.01 of the value at zero separation. It would be interesting to investigate these corre- 
lation functions in lattice calculations. In the derivation of the identity (11061) . we made use of 
the transformation properties of the field- strength under rotations by 90 degrees, and one may 
therefore expect that such relations remain valid on the lattice with the discrete symmetry of a 
hypercube. 

It will also be important subsequently to note the close relationship of T to the Poynting vector. 
Using the same arguments as above, it follows that the off-diagonal matrix elements T_ are 
related to the Poynting vector by the rotation (11031 ) 

P = E"xB°, 

with the consequence that 

3 

C T _(0) = (r_ y (0)r_y(0)> = £ e^(P fc (0)P z (0)) fixed. (107) 

k,l=l 

Wilson loop operators 

In addition to the local operators discussed above, it is also interesting to consider correlation 
functions of Wilson loops. This will enable us to investigate the role of the locality of the 
operator in the determination of the glueball masses and to make contact with the operators 
used in lattice gauge calculations. 

We denote the direction in which the Wilson loops are separated as time direction and will 
perform weighted spatial averages of the Wilson loops in order to project on states of definite 
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spin. We consider circular loops of radius r and denote the unit vector normal to the plane of 
the circle by n and the coordinate of the center of the circle by x , 

W r (n,xo,t) = IxrPe^^^^^. (108) 

If applied to the vacuum, this operator generates states of different momenta and different spins. 
Projection of W(n, x , t) on definite spin £, m is obtained by 

< m (x , t) = J dn n W r (n, x , t) F/ m (n) (109) 
and on zero momentum by spatial averaging 

W* m (t) = i J d 3 x ^ m (x ,t) . (110) 

Translational invariance in Euclidean 4-space and 3-dimensional rotational invariance impose 
constraints on the correlation functions. We have 

((W? m (t) - (W? m (t))(Wf m '(t') - (W* e ' m> '(*'))> 

= S ul 8 mml [{W l r {t - t')W?(0)) - (Wf) 2 6t, ] =C t (t- t') 5 u ,5 mm , . (Ill) 
We have used the fact that by rotational and translational invariance the correlator 

(/ d 3 x W,(n,x ,t) J d 3 y W r *(n',y Q ,t')) = W r (n-n',t-t') (112) 

can depend only on the angle between the two orientations of the Wilson loops. Similarly 

((^ m (x ,t) - (W^^ ,t)))(Wf m '( Xo ,t') - (W* e ' m '(x ,t))) = D e (t-t')6 u ,6 mm , . (113) 
We finally observe that [cf. Eq. (11081) 1 

W r (n,Xo,t) = W r (-n,Xo,t), (114) 
and therefore the correlation functions (II 13111 14) vanish for odd £. 

In order to improve the statistics, we actually do not make use of the m-independence of the 
correlation functions but rather sum over all values of m in Eqs. (|1 1 1|) and (II 12|) . In order 
to relate the Wilson loop correlation functions to the field- strength correlation functions, we 
consider the limit of small size Wilson loops, i.e., r p. In this limit and for fixed orientation 
n we have 

r^O: W r (n,Xo,t) = #Vr 4 (n ■ B a (x ,t)) 2 . (115) 
Integration over the orientation [Eq. (11091) 1 yields 



\/Air 

r - : ty r <°(x ,t) = ^W(B*(x ,t)) 2 , (116) 
whereas the £ = 2 projection gives, for instance, for the m = component 

r - : W r 2 '°(*o,t) = ^AV[(BJ(xb,f)) 2 - ^(B a (x , tf)} . (117) 
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As expected, at small sizes, the angular momentum projected operators reduce for t = to 
(B a ) 2 and for £ = 2 to B <g> B, 

(B ® B) tj = B?B°: - \-WB a 5ij = 1 -{T + - T.)y . (118) 

o Z 

Electric field operators do not appear since the Wilson loops considered are located in spatial 
planes (e 4 • A a = 0). 

For subsequent reference, all the operators used to calculate correlation functions, salient prop- 
erties, and the results of measurements are summarized in Table |5] 



state 


name 


operator 


SCj <-> X4 


figures (slope) 


0+ 


S 


1/4(E 2 + B 2 ) 




|2I](-4.0) 


0+ 


W™ 


/ dQ n W r {n,x ,t)Y* {n) 






1- 


v P 






■■■■ - 10) 


1+ 


v a 






[26](— 10.5 .... - 11) 


2+ 


W 2m 


J-dn n W r (n,x 0) t)y 2 * m (n) 




[27](— 7.5 .... -7.8) 


2+ 


T_ 


E <g E — B <g B 


P = E a x B a 


M(-),m\(-)M(-) 


2+ 


T + 


E <g E + B <g B 


T a 


[H(-8.0),|23](-) 


2+ 


T a 


E (g B + B (g E 


T + 


|2Q](-),[22](-8.3),|23](-) 


2" 




e abc [(B a x B b ) ® E c + E c ® (B a x B 6 )] 




|25](-9.3) 



Table 5: Summary of operators used to calculate glueball masses and the results of measure- 
ments. The slope, specifying the glueball mass, for each figure is shown in parentheses. 



5.2.3 Correlation Functions of Field-Strength Operators 

With the preceding definitions and discussion, we now present the results of our calculations of 
correlation functions in instanton ensembles. 

T+ 

We begin with the correlation function of the 2 + operator T + [cf. Eq. (| 1001) 1 in Fig.QjB For small 
separations, the shape of the correlation function (curvature) still reflects the size of a single 
instanton, even though the single instanton correlation functions vanishes. Now the instanton 
size appears at small separations via the interference of some "mean-field" with the gauge field 
of a single instanton. In the intermediate regime .2 < \t\ < .6 approximately, a common slope 
can be associated with the correlation functions in these various ensembles. In this particular 
case of the T + correlation function, we extract a slope of 8. For larger separations we observe 
strong fluctuations in the correlation functions. As is shown in the lower part of the figure, 
these fluctuations are connected to violations of the rotational symmetry. After reducing the 
diagonal components of the tensor by 1/3, all 6 components of the tensor are identical if the 
ensemble exhibits exact rotational symmetry. Comparison of the two figures strongly suggest 
that the differences in the correlation functions for the different ensembles for \t\ > .5 — .6 are 
a consequence of the violations of the rotational symmetry in the individual ensembles. 
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Figure 19: Summed correlation function of T + [Eq. (11001) 1 and its decomposition into its 6 
components (denoted by different symbols) T +i j, i < j in an ensemble with 50 instantons of 
variable size. 

As a further consistency check we also have compared correlation functions of ensembles which 
differ in the number of pseudoparticles. In Fig. [20] we have rescaled the values of t for the 



1 i 1 i 1 i 1 i i.' i '.i i 1 i 1 i 1 i r 
t > 




Figure 20: Demonstration of scaling of the T_ (red and black squares, upper curves) and T a 
(green and blue diamonds, lower curves) correlation function for ensembles with 50 and 200 
instantons of size p = 0.16. 

correlation function of the ensemble with 50 instantons by a factor 1 .44 as suggested by the 
analysis of the Wilson loops (cf. Table [3]). As is seen, the extracted slopes or masses scale as 
expected with yfa within 10% or better. 

S 

The intermediate states contributing to the correlation functions of the action density are + 
states. Difficulties arise in the numerical evaluation of this correlation function due to the large 
contribution from the intermediate ground state. The subtraction of the vacuum expectation 
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Figure 21: Action density correlation function C s (x) of ensembles with 50 instantons of size 
p = 0.16 (red squares), p = 0.04 (green diamonds), and 800 instantons (L = 4) of size p = 0.16 
(blue circles). The slope of the straight lines is —4.0. 

value [cf. Eq. (l88~l) 1 makes the correlation function C s particularly sensitive to the violations in 
translational and rotational symmetry. A stable and unambiguous result could be obtained only 
by constructing ensembles in a larger volume in which the violation of translational invariance 
is significantly weaker. We have used the ensemble of pseudoparticles discussed above in the 
context of the thermodynamic limit and calculated the action density correlation function in a 
16 times larger volume with 800 pseudoparticles. The result of this calculation is shown in 
Fig. El 




-1 -0,8 -0,6 -0,4 -0,2 0,2 0,4 0,6 0,8 1 



Figure 22: Correlation function of T a for different ensembles. 800 instantons (L = 4) of size 
p = 0.04 (blue diamonds), 50 instantons of size p = 0.04 (black squares) and of variable size 
(red circles). The value of the slope is —8.3. 

Unlike for the other correlation functions, very different results are obtained for the action den- 
sity correlation in the two "equivalent" ensembles. We note that the correlation function in the 
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small volume displays a zero at \t\ = 0.75 independent of the instanton size. In the larger vol- 
ume, these fluctuations violating translational invariance are reduced resulting in a reasonable 
behavior of the correlation function. As illustrated in Fig. [2T|for sufficiently small instanton 
sizes, where the single instanton contribution is restricted to small separations, a fiducial inter- 
val exists where a meaningful value of the mass can be extracted also in the ensemble with the 
smaller volume, and we indeed obtain the same slope. 

T+,T_,T a 

The correlation functions of the 3 tensor operators bilinear in the field-strength [cf. Eq. (| 1001) 1 




-0,8 -0,4 0,4 0,8 



Figure 23: Tensor correlation functions T + (blue circles), T_ (black squares), and T a (red dia- 
monds) of an ensemble with 50 instantons of size p = 0.1. 

with the 2 + (T±) and 2~ (T a ) excitations appearing as intermediate states are shown in Fig.l2~3l 
Surprisingly, not the correlation functions of operators with the same but rather those with oppo- 
site parity are very similar in shape and actually possess the same normalization. Our numerical 
results yield Ct + (0) = Cy a (0) within the statistical accuracy and furthermore, this equality per- 
sists at finite separation within a few percent for the corresponding values of the momentum 
projected (|99l ) correlation functions. As discussed in the previous section, the equality at zero 
separation follows from Eq. (11061 ), arising from the equivalence of the two operators under 
interchange of x t and Xj. 

We now turn to the extraction of glueball masses from these correlation functions and discuss 
first the case of the negative parity 2~ state excited by the tensor T a . As Fig. [22] shows, a 
sufficiently large fiducial volume exists for extracting masses for the two ensembles and the 
two values agree within 5% with an average of 8.3. We will come back to this state when 
discussing products of three field- strength operators. 

Similarly, the analysis of the correlation function of T + already discussed above (cf. Fig. IT9T ) 
yields for the same ensembles a common slope of about 8. The discussion of the Wilson loop 
correlation function will provide additional information on the corresponding 2 + state. 
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Also in the T__ correlation function, 2 + states appear as intermediate states. As the left side of 




Figure 24: Left: Correlation function of T_ for different ensembles (Nj, p) = (50, 0.16) - green, 
(50, 0.08) - red, (50, variable p) - indigo. Right: Correlation function of T_ (red) and compo- 
nents for the ensemble (50, 0.08). The values of the slopes are 4.5 and 6.8. 

Fig. [24] shows, this correlation function is rather exceptional, in that the fiducial intervals for the 
different ensembles extend to distances of up to 1. A fit to these correlation functions for \t\ < 
0.6, the fiducial interval also used for the other correlation functions, yields an average value of 
4.1 for the slope. At larger values of \t\, the slopes of the correlation functions for the different 
ensembles start to deviate. As is seen from the right side of this figure, the correlation function 
for the (50, 0.08) ensemble displays an increase in slope by a factor 1.5. Furthermore, as the 
decomposition in different components indicates, this change of slope is not accompanied by a 
significant violation of rotational invariance. A change to a larger slope also occurs though more 
moderately around t = 0.6 in the (50, 0.16) (from 3.7 to 4.9) and in the (50, variable p) (from 
4.1 to 5.2) ensembles. For determination of a glueball mass, we have to use the larger values 
of the slope as more realistic. The change to a larger slope with increasing separation suggests 
that for this particular operator, the regime in which the time evolution is governed by a transfer 
matrix has not been reached in the region \t\ < 0.6 where the smaller slopes are extracted. As 
for the T + - T a relation, our numerical results for the momentum projected correlation functions 
of P and T_ agree for \t\ < .4 within 6%, indicating that, like the correlation function of 
the integrated Poynting vector, the correlation function of the integrated T operator has no 
interpretation in terms of glueballs in this regime. Despite their rather complicated structure, the 
correlation functions of T_ exhibit the expected scaling behavior as shown in Fig.[20l However 
the strong dependence of the larger slope on the instanton size (4.9 - 6.8) prevents extraction of 
a glueball mass from this correlation function. Given the comparatively small fiducial interval 
(0.7 < \t\ < 1), we cannot rule out that several states contribute with weights depending on the 
instanton size. In this case we can only deduce an upper limit for the mass 

m 2+ <4.9. (119) 
We turn now to a discussion of correlation functions associated with observables composed of 
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products of 3 field- strength operators. We will study excitations of l 1 * 1 states and reexamine 
the excitation of 2 ± states. Since the field- strengths are evaluated at the same point, very rapid 
variations in the vicinity of the center of the pseudoparticles result. 




Figure 25: T 3 correlation function for the ensemble of 50 instantons of variable size (black 
circles) and 800 instantons (L = 4) of size p = 0.16 (red squares), the value of the slopes is 
-9.3. 

Due to these rapid variations significant results could not be obtained for the 0^ excitations 
which contain single instanton contributions. In the case of the 1* and 2 ± excitations, where 
single pseudoparticles do not contribute, the variations due to the interference with the "mean 
field" are still too large for generating reliable results for the smallest instanton size (p = 0.04) 
considered. Consequently the fiducial intervals are in general significantly smaller than in the 
case of products of two field- strength components. The correlation function of the 2 operator 
T 3 shown in Fig. [25] displays at small \t\ rapid variations and a strong dependence on the indi- 
vidual ensembles. In particular, the ensemble with variable p exhibits for small \t\ an almost 
singular behavior. Beyond this region, a sufficiently large fiducial interval exists where a com- 
mon slope can be extracted for the two ensembles. The value of 9.3 differs by 10% from the 
value 8.3 in Fig. [22} The average of these two values will be used subsequently. A meaningful 
extraction of a slope parameter of the correlation functions for the 2 + operator (E x E) <g> B 
was not possible and is not discussed here. 

V a ,V p 

We finally turn to a discussion of the spin 1 correlation function in Fig. [261 For both, the V a 
and V p [cf. Eq. (11021) 1 correlation functions, the three ensembles exhibit, beyond the small \t\ 
region, sufficiently large fiducial intervals for extracting the slopes reliably. The average values 
of these slopes are used later. 

For various ensembles, we have investigated Wilson loop correlation functions for loops of 
different radii. We have considered both Wilson loops of fixed radius and homogeneously dis- 
tributed radii. For an appropriate choice of the parameters, the observables have been found to 
be insensitive to this choice. Due to the higher computational complexity as compared to the 
correlation functions of the local operators, less configurations and for a given configuration 
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Figure 26: Axial V a (left) and polar vector V p (right) correlation functions for ensembles of 50 
instantons of variables size (black circles ), of size p = 0.16 (blues squares) and p = 0.08 (red 
stars). The values of the slopes are -9.0, -9.7, -10.0 (V a ) and -10.5, -11.0 (V p ). 

fewer points for the momentum projection had to be used. The use of Wilson loop correlation 
functions did not improve the determination of the + glueball mass, since also in this case dif- 
ficulties due to the subtraction of the vacuum expectation value with the concomitant sensitivity 
to deficits in rotational and translational invariance occur. As shown by the right side of Fig. [271 
slopes (or masses) could be extracted from the 2 + correlation functions with an average value 
of about 8 and an uncertainty of about 15%. Form this figure we also can read off the relative 
normalization of the correlation functions for the different sizes of the Wilson loops. As com- 
pared to the correlation functions for (r) = 0.2 the strength of the coupling to the 2 + excitations 
is enhanced by a factor of 5.7 for the correlation function with (r) = 0.4. One therefore might 
expect a larger systematic uncertainty in the slope extracted for (r) = 0.2. As the figure shows, 
the size of the Wilson loop that produces the maximal overlap with the exited 2 + state is 0.4 and 
the extracted slope essentially does not change when varying the size from 0.3 to 0.4. Adopting 
the rational that the maximal overlap produces the most meaningful slope, we extract from this 
correlation function the value 7.7. We observe that this result is within ~ 10% compatible with 
the slope extracted from the T + correlation function. We also note that the fact that the overlap 
is maximal for a Wilson loop of radius r = 0.4 implies that the size of the state is approximately 
0.47 fm, consistent with lattice calculations ||2~9ll3~0"l . 



5.3 The Glueball Spectrum 

The results for the glueball masses extracted from our investigations of correlation functions of 
local operators built from 2 and 3 field- strength operators and from the non-local Wilson loop 
operators are summarized in Table [6] After a rescaling, our results are compared in Fig. [281 to 
the results of the lattice gauge calculation of Ref. IT3TI . In comparison with the lattice results, 
the energy scale is too small by roughly a factor of 2. For a more detailed comparison we have 
rescaled our results with a factor of 2.2, which makes the value of the average excitation energy 
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Figure 27: Left: Wilson loop correlation functions D e (t) (11111) for £ = (red squares), £ = 2 
(blue triangles), iV = 50, p = 0.16 and (r) = 0.1 compared with field- strength correlators of B 2 
(black circles) and B ® B (green diamonds). The two £ = 2 correlation functions are virtually 
indistinguishable. Right: Wilson loop correlation functions D 2 (t) (II 111) for N = 50, p = 
0.08 and (r) = 0.2, 0.3, 0.4, 0.5 (green diamonds, blue circles, blacks stars, red squares). The 
extracted slopes are —6.6, —7.8, —7.5, —9.3. 



state 


slope 


slope/ y/a 


0+ 


4.0 


1.6 


2+ 


7.7 


3.1 


2~ 


8.8 


3.5 


1+ 


9.6 


3.9 


\- 


10.8 


4.4 



Table 6: Glueball masses. Note that errors are discussed in the text and are not indicated here. 



agree for the two spectra. Having adjusted the overall scale, we note that the masses of the 
+ , 2~ and 1 ± glueballs of the two spectra agree within 5% while the mass of the 2 + state is 
25% larger than the lattice value. We find this agreement surprising also in view of the large 
adjustment of the overall scale. These final results do not account for the possible existence of 
a second 2 + with an upper limit (after rescaling) for the mass of [cf. Eq. (11191) 1 

m 2 + < 2 y/a . 

As emphasized above, the smallness of the fiducial intervals did not permit us to firmly establish 
or rule out the existence of an additional low mass 2 + state. 

The correlation functions for meron ensembles yield a very similar spectrum although the un- 
certainties in the extraction of the slopes, for the reasons discussed above, are significantly 
higher. In particular, finer details like the splitting of the 1^ states could not be determined. 
As for the instanton ensembles, the existence of a 2 + state, essentially degenerate with the + 
excitation, cannot be ruled out. The overall scale factor necessary to reproduce the average 
excitation energy of the lattice results is 1.8. 
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Figure 28: Spectra of + ,2 ± , l± glueball states of SU(2) Yang-Mills Theory. Left: Lattice 
gauge theory [31]. Right: Masses from pseudoparticle ensembles after rescaling (see text). 



5.4 Wilson Loop Correlation Functions and Confinement 

In this concluding paragraph of our studies of correlation functions, we will return to the issue 
of confinement in the pseudoparticle ensembles. The physics of confinement is intimately re- 
lated to the formation of flux-tubes or gauge strings and one may wonder whether the area law 
found in the fit to the Wilson loops indeed implies that such a mechanism is operative also in 
the pseudoparticle ensembles. To study the dynamics of Wilson loops we follow the investi- 
gations in 11321 and calculate correlation function of loops with fixed orientation. We assume 
circular loops of equal radius which are parallel to each other and orthogonal to the direction 
of separation. We do not perform any averaging, i.e., we calculate the following correlators [cf. 
Eq.flDSfc] 

C r (t) = (W r (n, Xq, t)W r (n, x , 0)) . (120) 



The results of the calculations are shown in Fig. [291 If one assumes that the application of the 
Wilson loop operator to the vacuum generates a gauge string with (approximate) energy 



E = lixor 



the relevant variable for describing the correlation function is rt. As Fig. [29] shows, the three 
correlation functions are described approximatively by a universal curve in terms of this vari- 
able. It thus appears that the dynamics in the pseudoparticle ensembles is compatible with the 
formation of gauge strings or flux tubes. The existence of such a universal curve offers the pos- 
sibility for an alternative determination of the string tension. Due to the limited range of values 
a direct extraction of the string tension is not possible. We have parameterized the distributions 
by 

lnW = - ^ . (121) 

y/b + a(rt) 2 

From the fit to the correlation functions of the Wilson loops shown in the right part of Fig. 
we obtain 



-5^ = 17.4 ±6. (122) 
2ir 
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Figure 29: Left: Logarithm of correlation functions of circular Wilson loops with fixed parallel 
spatial orientation for an ensemble of 500 instantons of size p = 0.08 as a function of separation 
in time t for 3 values of the Wilson loop radius r = 0.25 (black circles), r = 0.32 (red squares), 
r = 0.40 (blue diamonds) . Right: The same as a function of the variable rt. 



The small range of relevant values of the variable rt does not permit a more precise determina- 
tion of a. The value (|122|) of the string tension agrees within 20% with that derived from the 
asymptotics of the Wilson loops (cf. Eqs. (|75T) and (|76|)). 

At large separations the universality breaks down as expected. For rt > r 2 it becomes energet- 
ically favorable for the system of two loops to couple by annihilation of the gauge string to the 
vacuum or glueball states. This corresponds to the Gross-Ooguri phase transition in the large 
N c limit E3. 



6 Conclusions 

The effective theories with pseudoparticle degrees of freedom studied in this work describe im- 
portant properties of QCD. The regular gauge instantons or the merons with their long-range 
gauge fields apparently constitute a proper and economical set of degrees of freedom in the de- 
scription of the Yang-Mills dynamics in the infrared. By superposition of these building blocks, 
field configurations are generated with confinement as an inherent property. An area law for suf- 
ficiently large Wilson loops is already obtained in (stochastic) ensembles with the same weight 
assigned to all field configurations. However such ensembles violate basic principles. The 
action of field configurations in such ensembles is diverging in the infrared, concomitantly the 
ensembles do not possess a proper thermodynamic limit and violate translational invariance. By 
assigning the standard weight given by the action to the field configurations, all these problems 
are cured at once and confinement is preserved. The resulting configurations not only confine 
static color charges but also the building blocks themselves. Only at the expense of a logarith- 
mically diverging action can a building block be removed to infinity. The pseudoparticles are 
strongly correlated over large distances in order to prevent a logarithmic rise in the action (for 
this reason a global update procedure had to be employed in our computations). On the other 
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hand with its many spikes and valleys in color singlet quantities like the action density or the 
topological susceptibility, single field configurations appear strongly disordered. 

Thus, in our investigations, a picture of the Yang-Mills vacuum emerges which is that of a 
nematic substance - a medium that is strongly correlated in color or "director" space (S 3 and 
RP 2 respectively) and which at the same time is essentially disordered in configuration space 
(cf. lT34lO . This picture summarizes our findings. On the one hand, updates of the randomly 
chosen coordinates of the pseudoparticles in addition to the color update have no significant 
effect. On the other hand, a color disordered ensemble such as the stochastic ensemble we 
have studied does not produce a qualitatively correct model of the Yang-Mills vacuum. Clearly, 
the color correlations in meron and instanton ensembles are at the heart of the phenomenon of 
confinement. When a singular gauge instanton contributes to 2 Polyakov lines, it produces a 
linear potential, while singular instantons affecting one Polyakov line only are irrelevant. The 
last part of this argument [5 ] gets modified in our context. Due to the strong color correlations, 
single instantons still contribute to the potential. In addition to confinement of static charges, 
the color correlated fields reproduce remarkably well other important features of the Yang-Mills 
dynamics. They exhibit proper scaling under variations of the number of pseudoparticles or the 
size of the system and they reproduce qualitatively characteristic properties such as the values 
of the gluon condensate and the topological susceptibility. Unlike the gluon condensate, the 
value of the topological susceptibility can be related to the single pseudoparticle properties. 
We also have been able to investigate finer details and for instance to show analytically that 
characteristic properties of the Wilson loops in the limits of small and large size, i.e., the U(l) 
and the string limit, are qualitatively reproduced. Here the role of the pseudoparticle size as 
an effective ultraviolet regulator is essential. We also qualitatively confirmed the string picture 
for Wilson loop correlation functions and thus obtained a significant consistency check of the 
description of confinement in term of pseudoparticles. 

The outcome of our investigations of correlation functions is more ambiguous. We have demon- 
strated the emergence of a hadronic scale under the usual circumstances where the size of the 
building blocks is smaller than the emerging size of the "many-body" hadron. In the case of 
merons with their long range field-strength the opposite mechanism is at work and establishes 
a hadronic size smaller than the constituent meron size. To study hadronic excitations, we have 
computed momentum projected correlators of local gauge invariant operators with definite an- 
gular momentum built from field-strength bilinears and trilinears and have also calculated cor- 
relation functions of momentum and angular momentum projected Wilson loops. With these 
computations, we have been able to determine the spectrum of glueballs with quantum numbers 
+ ,l ± and c l ± . While the structure of the spectrum agrees qualitatively with the correspond- 
ing lattice results, the masses determined are too small by an overall factor of about 2 for both 
instanton and meron ensembles. We have not been able to uniquely identify the origin of this 
difficulty. Ambiguities in the determination of the string constant are unlikely to explain it. A 
smaller value deduced from the fit to the Wilson loop would, after identification with the em- 
pirical value, result in a higher energy scale. There is indeed some freedom left by redefinition 
of the relevant interval of the area of the loops. A change of the order of 20% cannot be ruled 
out entailing a 10% change in the energy scale. Corrections of a similar order of magnitude 
seem to result when using a "smeared" Wilson loop [35]. Neither effect however approaches 
the necessary factor of 4 in the string tension. The most likely explanation is that with the pseu- 
doparticle ensembles primarily determined by infrared properties, the ultraviolet fluctuations 
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are not adequately described. For a description of an object like the scalar glueball with a size 
of approximately 0.2 fm [29] the spatial resolution of our field configurations - possibly limited 
by a too "small" number of pseudoparticles - may just be insufficient. In our discussion of small 
size Wilson loops we indeed have noticed the disappearance of the "Coulomb regime" with the 
disappearance of the disorder in color space. To cure this problem one would have to either sig- 
nificantly enlarge the number of pseudoparticles or to introduce additional degrees of freedom 
which improve the description of the vacuum at small scales. Candidates for such additional 
degrees of freedom are singular gauge instantons which are the center reflected partners of the 
regular gauge instantons. Inclusion of these degrees of freedom would, by construction, avoid 
the otherwise problematic issue of multiple counting of gauge copies. 

In addition to such improvements in the description of the Yang-Mills theory, for phenomeno- 
logical applications quarks have to be incorporated. Here the advantages of the continuum for- 
mulation of the effective pseudoparticle theories could be essential. First steps in this endeavor 
have been taken [36]. Here one aims to take over the successful treatment of the dynamics of 
quarks in the instanton liquid model [H. The ultimate goal would be to establish in this way an 
effective theory that exhibits the two nonperturbative basic phenomena of QCD - confinement 
and chiral symmetry breaking. 
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